Hecke Algebras, SVD, and Other Computational Examples 

with CLIFFORD* 

^— > . 

Rafal Ablamowicz 
>-' ■ Department of Mathematics, 

Tennessee Technological University 
Cookeville, TN 38505 
"""J 1, rablamowicz@tntech.edu 

<: 

Ql ■ September 30, 1999 

I: 

S Abstract: CLIFFORD is a Maple package for computations in Clifford algebras C£(B) of an arbi- 

trary symbolic or numeric bilinear form B. In particular, B may have a non-trivial antisymmetric 
part. It is well known that the symmetric part g of B determines a unique (up to an isomorphism) 

t-H ■ Clifford structure on C£{B) while the antisymmetric part of B changes the multilinear structure of 

C£(B). As an example, we verify Helmstetter's formula which relates Clifford product in C£(g) to 
the Clifford product in C£(B). Experimentation with Clifford algebras C£(B) of a general form B is 
highly desirable for physical reasons and can be easily done with CLIFFORD. One such application 

^^ ■ includes a derivation of a representation of Hecke algebras in ideals generated by q -Young operators. 

t-H , Any element (multivector) of C£(B) is represented in Maple as a multivariate Clifford polynomial in 

the Grassmann basis monomials although other bases, such as the Clifford basis, may also be used. 
Using the well-known isomorphism between simple Clifford algebras C£(Q) of a quadratic form Q 

,_£^ ' and matrix algebras through a faithful spinor representation, one can translate standard matrix al- 

gebra problems into the Clifford algebra language. We show how the Singular Value Decomposition 
of a matrix can be performed in a Clifford algebra. Clifford algebras of a degenerate quadratic form 
provide a convenient tool with which to study groups of rigid motions in robotics. With the help 
from CLIFFORD we can actually describe all elements of Pin(3) and Spin(3). Rotations in R 3 
can then be generated by unit quaternions realized as even elements in C£q 3 . Throughout this work 
all symbolic computations are performed with CLIFFORD and its extensions. 
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1 Introduction 

A first working version of a Maple package CLIFFORD was presented in Banff in 1995 (D. From 
a modest program capable of symbolic computations in Clifford algebras of an arbitrary bilinear 
form, CLIFFORD has grown to include 96 main procedures, 21 new Maple types, close to 4,000 
lines of code written in the Maple programming language, and an extensive on-line documentation. 
There is a number of special-purpose extensions available to CLIFFORD such as SUPPL and ASVD 
used in this paper [n. In fact, anyone who uses Maple can easily write additional procedures to 
tackle specific problems. 

There are major advantages in using CLIFFORD on a Computer Algebra System. One is an 
ability to solve equations and find the most general elements in the Clifford algebra satisfying given 
conditions. This approach has been presented in Sections and H. In Section Young operators 
in the Hecke algebra Hf(3,q) are eventually found by systematically solving three equations that 
define them. Computations in this section were first reported in W along with a physical motivation. 



There, experimentation with CLIFFORD led to finding Young operators realized as idempotents in 
the Hecke algebra which in turn had been embedded into the even part of a Clifford algebra C£(B) 
of a suitable bilinear form B. Such embedding was first given in ]17| where the bilinear form B was 
found so that the defining relations on Hecke generators were satisfied. Furthermore, it was shown 
in Q that Young operators corresponding to conjugate tableaux were related through the operation 
of reversion in the Clifford algebra, an idea that was first proposed in flTJ]. Here, we only show the 
mechanics of the search for such operators, Garnir elements, and bases in the representation spaces 
as they were performed with CLIFFORD. 

In the same spirit, in Section g, we describe a search for the elements in Pin(3) considered as 
a subgroup of the group of units in the Clifford algebra C^o,3- Seven general types, not entirely 
exclusive, are eventually found through a systematic search and analysis. Then, the elements of 
Spin(3) are computed and related to unit quaternions. Rotations in coordinate planes and in a 
plane orthogonal to an arbitrary non-zero axis vector are described using quaternions realized as 
elements of C£q 3 . A symbolic formula describing the most general rotation is derived. Finally, 
using the ability of CLIFFORD to compute in Clifford algebras of a degenerate quadratic form, 
the semi-direct product Spin(3) x M 3 is shown to generate rigid motions on a suitable subspace of 
the Clifford algebra C£q,3,i- 

The second advantage of using symbolic program like CLIFFORD is its ability to compute 
with expressions containing totally undefined symbolic coefficients. It is possible, of course, like 
in Section H to impose additional conditions on these coefficients when needed (by defining aliases 
for roots of polynomial equations). In Section H we verify one of Helmstetter's formulas pj| that 
relates Clifford product in C£(B) , the Clifford algebra of an arbitrary bilinear form B, to the Clifford 
product in C£{g) where g — g T is the symmetric part of B. A re- wording of the Helmstcttcr formula 
presented to the Author by Pertti Lounesto pj| proved to be suitable for symbolic verification with 
CLIFFORD. In fact, while this problem turned up to be a challenge for CLIFFORD in view of 
its complexity, it also has helped to fine-tune the program to make such computations feasible. We 
will only illustrate computations in dimension 3; however, computations in dimension up to 9 have 
been successfully completed. 

The third major advantage of using CLIFFORD shows up in Section H: here, we perform within 
the same workspace not only symbolic computations with a Clifford algebra, but also with a linear 
algebra package built into Maple. This way we can illustrate in two low-dimensional examples two 
parallel approaches to the Singular Value Decomposition (SVD) of a matrix: one through the matrix 
algebra, and one through the Clifford algebra. In Section [h| we comment on SVD based entirely 
on the Clifford algebra approach. 

2 Verification of the Helmstetter formula 

In his paper |jl| Helmstetter studies canonical isomorphisms between Clifford algebras C£(Q) and 
C£(Q') of two quadratic forms Q and Q' defined on the same (real or complex) vector space V. 
The forms are related via the identity Q'(x) = (J(x) + B(x, x) for every x £ V and some bilinear 
form B on V. Helmstetter constructs a deformed Clifford product * on C£(Q) by extending the 
Clifford product xy of two elements x and y in V '—* C£(Q) 

x*y = xy + 5(x,y) 



to all elements in Cliff). Together with the new product *, the Clifford algebra C£(Q) becomes 
a deformed Clifford algebra C£(Q,B). Given now two different bilinear forms B and B' on the 
quadratic space (V, Q) such that ~B(x, x) = -B'(x, x) for every x G V, Helmstetter proves that there 
exists F € /\ V such that 

B / (x,y)-B(x,y)=<F,xAy> 

and the mapping 

4>:C£(Q,B)->C£(Q,B'), u ^ e AF J u (1) 

gives an isomorphism from C£(Q, B) to C£(Q, B') which acts as an identity on V. In the above, 

e AF J u denotes the left contraction of u by the exterior exponential of F (see ||, J22|). A special 

case of (Q) occurs when B is symmetric, that is, B = g = g T and B' = g+A for some antisymmetric 

form A. 

Thus, with a slight change of notation, let B = g + A, g T = g, A T = —A and let us consider 

two Clifford algebras C£(g) and C£(B) on the same vector space V. We have therefore three 

contractions: x J y = B(x, y), x J y = j4(x,y), and x J y = <?(x, y).p] Then, the B -dependent 
B A ' g "- 1 

Clifford product uv of any two elements u and v in C£{B) can be written |23J ] solely in terms of 

B 

the operations in C£(g) as 

TO = ((« L e AF )(w L e AF )) L e A (- F ) (2) 

B g g g 

where w L e AF denotes the right contraction of u by e AF with respect to g. The product of 
g 

u L e AF and v L e AF in (0) is taken in C£(g). Let dimR(V) = n. Then, the element F g A V is 

g g 

dehned as 

F = J2(-ir iK - L)l A K e L r 1 (3) 

K,L 

where j _1 denotes the inverse of the unit pseudoscalar j = ei A e2 A ■ • • A e„ in C£(g), the product 
etj -1 is taken in C£(g), and the summation is taken over all multi indices K = [kijkz] and 
L = [h, I2, ■ ■ ■ , l s ] satisfying the following relations: 

K n L = 0, K U L = {1, 2, . . . , n], n = 2 + s, K and L are ordered by < . 

tt(K,L) denotes a permutation which puts the list [ki, k%, h, I2, ■ ■ -,1s] m the standard order 
[l,2,...,n] and \n(K, L)\ equals or 1 depending whether ir(K,L) is an even or odd element 
of S n . In (g) we have also adopted notation e^ = ej 1 ; 2 ...; s = e; A ei 2 A •• • A e; 3 . 

Before we proceed to verify formula (g) with CLIFFORD B, let's observe the following prop- 
erties of the left and right contraction: 

ut v =r 1 (du)Av), u J v = (uA(vj))j- 1 , J _1 = TTTT' ( 4 ) 

g g det{g) 



1 The symbols x J y, x J y, and x J y denote the left contraction of y by x with respect to B, A, and g 
B A g 

respectively. 



where ~ is the g -dependent reversion in C^(g).n Observe also that since F <E /\ V, we have 



oA-F 



7>Afc 



N T?r 
S— ' iV=L " /2j 

k=Q 



(5) 



where F Afc — F A F /\ ■ ■ ■ /\ F is the exterior product of F computed fc -times and ['J denotes the 
floor function. For example, for different values of n, F has the following form: 



n = 2, F = -Ai 2 e 12 1 , 

n = 3, F = -(A u e 3 - A i3 e 2 + A 23 ei) e'g, 
n = 4, F = -(A12 e 34 - Ai 3 e 24 + ^4i4 e 23 + A 23 e i4 + ,4 34 e i2 - ,4 24 e i3 ) e^ 



1234^ 



(6) 

(7) 
(8) 



and so on. We will verify the validity of i 
language, formula (0) becomesfl 



in a numeric and a symbolic case. In the Maple symbolic 



cmul(u, v) = RCg(cmulg(RCg(u, wexp(F, iV)), RCg(w, wexp(F, N)), wexp(-F, N))) 

with the CLIFFORD procedures cmulg and RCg representing the Clifford product and the right 
contraction in C£(g) and wexp giving the exterior exponential in /\ V. We limit our two examples 
to dimension n — 3. Computations presented in the following two sections can be extended with 
CLIFFORD to higher dimensions. 



2.1 Numeric example when n = 3 

Let's first assign an arbitrary matrix to B, split B into its symmetric and antisymmetric parts g 
and A, and compute the bivector F with a procedure makeF: 

> dim: =3: eval(makealiases(dim, 'ordered' ) ) :B :=matrix(dim,dim, [4,8,3,0,9,5,-2,1,7]) ; 



> g,A:=splitB(B) ; 



9, A:-- 



B : = 


l 

( 


1 8 


3 

5 






-2 1 7 




[4 4 il 
2 







4 * 
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4 9 3 
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-4 
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U' 


7 




-5 

- ~Y 


-2 



> F : =makeF (dim) ; 



8(» 



455 



el2 



F := —el3 
91 
Next, we find the exterior exponentials of F and — F which we assign to Maple variables Ft and 

F 2 respectively. 

> N:=floor(dim/2) :Fl:=wexp(F,N) ;F2:=wexp(-F,N) ; 

2 From now on in this section we assume that det(g) ^ 0. 

3 In CLIFFORD, the Clifford product uv of two elements u and v can be entered as u &c v (the infix form) or 
as cmul(u, v) . 



2 8fi 2 Hfi 

Fi := 1 + — el3 + — el2, F2 := 1 - — el3 - — el2. 
91 455 91 455 

Let u and v be two arbitrary elements in C£(B) : 

> u : =2+el-e23+el23 ; v : =3-e3+el2+e23 ; 

u ■- 2 + el - e23 + el 23, v := 3 - e3 + el 2 + e23. 
The Clifford product uv of u and v in C£(B), the left hand side of (||), is then equal to 

B 

> cmul(u, v) ; 

-48 e3 + 79 Id + 13 el3 - 6 el2 + 81 e2 - 8 ei23 - 81 ei 
while the right hand side of (||) gives the same result: 

> RCg(cnrnlg(RCg(u,Fl) ,RCg(v,Fl) ) ,F2) ; 

-48 e3 + 79 Id + 13 e!3 - 6 el2 + 81 eg - 8 ei23 - 81 ei 



2.2 Symbolic computations when n = 3 

A purely symbolic computation when n = 3 will look as follows. Matrix B is now defined as 

an arbitrary symbolic 3x3 matrix with a symmetric part g and an antisymmetric part A, and 

F is again computed using the procedure makeF. The exterior exponentials of F and — F are 

again denoted respectively by F\ and F% . All symbolic parameters in B are assumed to be real or 

complex. 

> dim: =3: eval(makealiases(dim) ) : 

B : =matrix (dim , dim , [gl 1 , gl2+A12 , gl3+A13 , gl2-A12 , g22 , g23+A23 , gl3-A13 , g23-A23 , g33] ) ; 

gll gl2 + A!2 gl3 + A13 

13 : gl2 - A12 g22 g23 + A23 

glS - A13 g23 - A23 g33 



g,A:=splitB(B); 










.9, A := 
F: =map (normal , clicolle 


'gll gl2 gl3~ 
gl2 g22 g23 
gl3 g23 g33 

;ct (makeF (dim) ) ) 


5 
> 



-A12 
-A13 


A12 A13 
A23 
-A23 



(A23 gll + A12 glS - A13 g!2) e23 (-A23 g!3 + A13 g23 - A12 g33) el 2 



%l 

(A23 gl2 - A13 g22 + A12 g23) e!3 
+ %1 

%1 := -g33 g22 gll + g22 g!3 2 - 2 g!3 g23 gl2 + g33 gl2 2 



%l 



923' gll 

> N :=floor (dim/2) :F1 :=wexp(F,N) :F2 :=wexp(-F,N) : 

We will now define two general elements u and v in C£(B) by decomposing them over a Grassmann 
basis (provided by a procedure cbasis). Coefficients in these two expansions are assumed to be real 
or complex.^ 

> cbasis (dim) ; 



4 For technical reasons, in Maple these coefficients cannot be defined as ui and Vi\ that is why we use here um 
and Wi. 



[Id, el, e2, e3, e!2, el3, e23, el23] 

> u : =add (uu [k] *cbasis(dim) [k] ,k=l . .2 "dim) ; 

u := uui Id + UU2 el + uu^ e2 + uu^ e3 + UU5 el2 + uuq elS + uu-j e23 + uug, e!23 

> v:=add(vv [k] *cbasis(dim) [k] ,k=l . . 2~dim) ; 

v := vvi Id + W2 el + OT3 e2 + W4 e3 + W5 el2 + vvq el3 + vv? e23 + vvg el23 
The Clifford product of u and v in C£(B) is then collected and assigned to a constant resi which 
we won't display due to its length. 

> resl :=clicollect (cmul(u, v) ) : 

As before, we finish by computing the right hand side of (|2|). By assigning it to res-i, we can then 
easily find that res\ — resi = as expected. 

> res2:=clicollect(RCg(cmulg(RCg(u,Fl) ,RCg(v,Fl)) ,F2)) :resl-res2; 



3 Hecke algebra computations 

In |7J it was shown that the symmetric group S n and its group deformation, the Hecke algebra 
Hw(n,q) could be constructed as a subalgebra of a Clifford algebra C£(B) for a suitably cho- 
sen q -dependent non-symmetric bilinear form B. q -Young operators were constructed as Clifford 
idempotents and the Hecke algebra representations in ideals generated by these idempotents were 
computed. Appropriate q -Young diagrams and tableaux representing symmetrizers, antisymmetriz- 
ers, and operators of mixed symmetries were realized inside the Hecke algebra, while the ordinary 
case of the symmetric group was obtained in the limit q — ► 1. 

The Hecke algebra is the generalization of the group algebra of the symmetric group S n by 
adding the requirement that transpositions of adjacent elements i,i+ 1 are no longer involutions. 
Following M we set tf = (1 — q)ti + q which reduces to sf = 1 in the limit q — > 1. The defining 
relations of the Hecke algebra will be given according to Bourbaki JlM. Let {1, ti, ■ ■ . , t n } be a set 
of generators which fulfill these relations: 

ti = 0—q)ti + q, (9) 

titj = tjti, \i-j\>2, (10) 

tiU+lti = ti+lUti+l, (11) 

then their algebraic span is the Hecke algebra Hp(n,q). The algebra morphism p which maps the 
Hecke algebra into the even part of an appropriate Clifford algebra was found in |17[| . In particular, 
p(U) = bi := ej A ej+ n , i = 1, . . . , n, where ei, . . . , e2„ are the generators of the Clifford algebra 
C£(B,V), V = spanjei}, with the following non-symmetric bilinear form B : 

> dim : =8 : n : =dim/2 : eval (makealiases (dim , ' ordered ' ) ) : B : =def B (dim) ; 
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Then, the form of B guarantees that the following relations hold: 



b\ {l-q)bi + q, 

bjbi, \i-j\>2, 
bib l+1 bi = b i+ ibib l+1 . 



bibj 



(12) 
(13) 
(14) 



Following p7| , we define the Hecke generators bi as balanced basis Grassmann monomials of or- 
der 2 : 

> for i from 1 to n do b.i:=(e.i) few (e.(n+i)) od: 

Using procedure cliexpand we can expand these generators in terms of the unevaluated Clifford 
product which is denoted in Maple as &C : 

> seq(cat ( 'b' , i)=cliexpand(b. i) , i=l . .4) ; 

bl = {elkCe5)-qId, b2 = (e2 &C e6) - qld, b3 = (eSfeC e7) - qld, bi = (e4 &C e8) - qld. 
Checking if the Hecke generators satisfy the defining relations can be done as follows: 

> map(evalb, [seq(simplif y(CS(cmul(b. i,b. i)-(l-q)*b. i-q) )=0,i=l . .n)] ) ; 

[true, true, true, true] 

> map(evalb, [seq(seq(CS(cmul(b. i ,b. j )-cmul(b. j ,b. i) )=0, i=j+2 . . n) , j=l . . n)] ) ; 

[true, true, true] 

> map(evalb@f actorSnormal, [seq(CS(cmul(cmul(b. i,cat(b,i+D) ,b.i)- 
cmul(cmul(cat(b,i+l) ,b.i) ,cat (b, i+1) ) )=0, i=l . .n-1)]) ; 

[true, true, true] 
Let's define the remaining basis elements of the Hecke image p(Hp(4, q)) in CUE) : 



=CS(cmul(b2,b3)): 
=CS(cmul(b4,b3)): 
b234:=CS(cmul(b23,b4)) : 
b4321:=CS(cmul(b432,bl)): 



> bl2:=CS(cmul(bl,b2)) :b21 :=CS(cmul(b2,bl) ) :b23: 
b32:=CS(cmul(b3,b2)) :b34:=CS(cmul(b3,b4) ) :b43: 
bl21:=CS(cmul(bl,b21)) :bl23 :=CS(cmul(bl2,b3) ) 
b321:=CS(cmul(b32,bl)) :b432 :=CS(cmul(b43,b2) ) 
bl234:=CS(cmul(bl23,b4)) : 

Thus, p defines a homomorphism from the Hecke algebra H^(A, q) into the Clifford algebra C£(B). 
It was shown in |Tfl that p is not injective for n > 4, and that its kernel contains all Young diagrams 
which are not L -shaped. 



3.1 Hecke Algebra H v ( 2, q) 

We begin with the Hecke algebra Hp(2, q) generated by {Id, b{\ which reduces to S% in the limit 
q — > 1. We have thus only one q -transposition b\ from which we can calculate a g-symmetrizer 



i?(12) and a q -antisymmetrizer C(12). One of the features of the construction presented in M is 
that -R(12) and C(12) are related by the reversion in C£(B) : 

> R.12:=CS(q*Id+bl) ; #S_2 symmetrizer 
C.12:=CS(Id-bl) ; #S_2 antisymmetrizer 
evalb(reversion(C.12)=R.12) ; 

R12 := el5 + qld, cl2 := Id - el5 

true 

Observe, that -R(12) is the reversion of C(12), that is, i?(12) = C(12)~ where ~ is the reversion in 

the Clifford algebra Ct\ f \. Operators C(12) and i?(12) are almost idempotent and they annihilate 

each other: 

> factor (CS(R. 12 &c R.12)); 





(q + 1) {q Id + e!5 


(CSCC.12 &c C.12)); 






(q + l)(Id-e!5) 


c C.12, C.12 &c R.12; 





> R.12 

0, 
Upon normalization, the symmetrizer -R(12) becomes the q -Young operator Y^IJ while the anti 

symmetrizer C(12) becomes Y t 2 

mutually annihilating idempotents adding up to the identity element 

> Y2.12 := R.12/(l+q);Y11.12:= C.12/(l+q); 

qld+e!5 .. Id - e!5 



K2) 
,2 

The following computation verifies that Y± 2 and Y 1 2 are 



Y212 



Y1112 



q+l q+1 

ybas:=[Y2.12,Y11.12] :f :=(i, j)->CS(cmul(ybas [i] ,ybas[j])); 
YM:=evalm(linalg [matrix] (2,2, (i, j)->f (i, j) ) ) ;CS(Y11 . 12+Y2. 12) 

/ := [h i) -* CS(climul(2/&as 4 , ybasj)) 

'Id el5 

+ 



YM 



1 Q 





1 







Id 



e!5 



1 q+l 



Id 



In our construction, the Hecke algebra -£/f(2, q) is a subalgebra of the even part Ct[ x of Gl\^\ 
and it is generated by {ei,e5J. That is, the bilinear form G on the vector space spanned by 
{e 1; e 5 } is 

> G:=matrix(2,2, [B[l , 1] ,B [1 ,5] ,B [5, 1] ,B [5,5] ] ) ; 

( 

1 ( 

After symmetrization, diagonalization, and the limit q — > 1, G becomes diag(l, —1). Therefore, due 
to the isomorphism C£(B) ~ C£(g) (as associative algebras) we can view Hr(2, q) as the subalgebra 
of Gtii. This embedding implies that any idempotent of the Hecke algebra -£/f(2, q) must be an 
even Clifford element. The following computation shows that the Young operators found above are 
the only two nontrivial mutually annihilating idempotents in H^(2,q). 



G:= 



> xx:=CS(a*Id+b*bl) ; 

xx := a Id + b el5 

> sol : =clisolve2(CS(cmul(xx,xx)-xx) , [a,b] ) ; 

sol := [b = 0, a = 0, a = 1, b = 0, b = , a = , b = , a = ] 

g + 1 q + 1 g + l g + l 

> f . 1 :=normal (subs (sol [4] ,xx) ) ; f .2 :=normal (subs (sol [3] ,xx) ) ; 

qld+el5 M Id-elS 
g+ 1 g+ 1 

In this case, the Young operators happen to be the two even primitive idempotents fa, fa in Cl\ x . 
Notice, that in the case when q = 1, the idempotents /i and fa reduce to the well-known primitive 
idempotents 

> fll:=subs(q=l,fl); 

> f22:=subs(q=l,f2); 
in the Clifford algebra Ct\\. 



ill := -Id + -el5 

J 2 2 



f22 ---Id-- el5 

2 2 



3.2 Hecke algebra //f(3, q) 

In |Q Young operators related to various symmetries were constructed as idempotent elements in 
the Hecke algebra H^(3, q) embedded into the even subalgebra CK|" 2 of Ct^fl- As one of the main 
features of this construction, the Young operators corresponding to conjugate Young tableaux in the 
sense of McDonald |24j were related through the reversion in the Clifford algebra. The goal was to 
find four Young operators known to exist from the general theory of the Hecke algebras for n = 3. 
The four g -Young operators found had one parameter and generalized the four Young operators of 
S3 described in Hamermesh p| on p. 245. One of them was a full symmetrizer, one was a full 
antisymmetrizer and two were of mixed symmetry. 

The construction began with finding the most general element X in H^(3,q) such thatFl: 

X + X = Id (15) 

X 2 = X (16) 

XX = (17) 

In the first step, the most general element in the Hecke algebra H$(3,q) was found that satis- 

fled ©ft 

> bset: = [Id,bl,b2,bl2,b21,bl21] :X:=bexpand(add(K[i] *bset [i] ,i=l. .6)) : 

sollist[l] :=clisolve2(X+reversion(X)-Id, [seq(K[i] ,i=l. .6)] ) : 

X: =bexpand(subs(sollist [1] , X) ) ; 

5 See formulas (24), (25) and (26) in 0. 

6 Procedure clisolve2 is capable of solving the equation X + X — 1 = for the free parameters K\, . . . ,Ke 
appearing in X. However, when used repeatedly, it randomly selects free parameters in the solution. 
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X : = (-i K 6 q + X - K 3 q + 1 K 6 q 2 + 1 K 2 q + 1 - 1 A' 3 - 1 JsT 2 ) Id + A 2 M + X 3 ^ + #4 M2 

+ (-K 6 +K 6 q- K 4 ) b21 + K 6 bl21 
In the above, procedure bexpand expands elements in the Hecke algebra, which are normally ex- 
pressed in the Grassmann basis of Cl.2,2, in terms of the Hecke basis {1, &i, b 2 , &12, &21, &m}- 

Thus, the solution to ([15]) gives an element that belongs to a family parameterized by four real or 
complex parameters. In order to simplify Maple output, we define two aliases and then we substitute 
X found above into the second equation (O) . 

> alias (alpha=RootOf ( (1+q) *_Z~2+ (-q~2*K [4] +K [4] +q*K [2] -1+K [2] ) *_Z+K [4] *K [2] + 
K [2] ~2-K [4] -K [2] +q*K [4] -q~2*K [4] ~2-q*K [4] ~2-q~2*K [2] *K [4] +q*K [2] ~2) ) : 

> alias (kappa=RootOf ( (1+q) *_Z~2+ (-q~2*K [4] +K [4] +q*K [2] +1+K [2] ) *_Z+K [4] * 

K [2] +K [2] "2+K [4] +K [2] -q*K [4] -q~2*K [4] ~2-q*K [4] ~2-q~2*K [2] *K [4] +q*K [2] ~2) ) ; 

/, a, k 

> sollistxxl:=clisolve2(cmul(X,X)-X, [seq(K[i] ,i=l. .6)]) : 

xxlist := [seq (bexpand (subs (sollistxxl [i] ,X) ) , i=l . .nops(sollistxxl) )] : 
for i from 1 to nops (xxlist) do r . i :=bexpand (xxlist [i] ) od: 

The six representatives r\ , . . . , r$ defined in this last Maple command are the same as the elements 
displayed in nj\ after formula (30). Using the command f indbasis we can determine that the set 
{rj}, i — I, ... ,6, is of rank 4 : 

> nops(f indbasis ( [r . (1 . .6)] )) ; 

4 
Since the representatives are related by the reversion, we seek four linearly independent elements. For 
example, r%, r%, r 3 and r§ are linearly independent. From now on we assign them to /i,/2j/3,/4- 

> nops (f indbasis ( [xxlist [1] , xxlist [2] , xxlist [3] , xxlist [5] ] ) ) ; 

4 

> f 1 : =r . 1 : f 2 : =r . 2 : f 3 : =r . 3 : f 4 : =r . 5 : 

We can easily verify that the four elements {/1, f 2 , / 3 , /*} satisfy equation (Ua): 

> CS(f .l+reversion(f . 1)-Id) ,CS(f .2+reversion(f .2) -Id) , 
CS (f . 3+reversion(f . 3) -Id) , CS (f . 4+reversion(f . 4) -Id) ; 

0, 0, 0, 

and equation (|l6|): 

> simplify (CS(f.l &c f.l - f . 1) ), simplify (CS(f .2 &c f .2 - f .2)), 
simplify (CS(f .3 &c f.3 - f . 3) ), simplify (CS(f .4 &c f .4 - f .4)); 

0, 0, 0, 
Each of the four non-primitive idempotents fi generates a three-dimensional one-sided ideal in 
H$(3,q). In order to split such ideal into a one-dimensional space and a two-dimensional space, 
one had to find a way of splitting at least one of these idempotents into a sum of two mutually 

annihilating idempotents. It was observed in pj that /1 contained a full symmetrizer. Since the 

(3) 
full symmetrizer Y 1 2 3 was defined there as 



,(3) _ g 3 l + g 2 bi + q 2 b 2 + gb 12 + gb 2 i + 6121 
< 2 ' 3 : ~ (l + q + q 2 )(l + q) 



v {3) _ H J- ~r f <n-TH "2 T <j"l2 ~r qu 2 i T 0121 ,-„■> 
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the full antisymmetrizer 7 123 was, by construction, the reversion of the symmetrizer, that is, 

2 1,2,3 ' — J l,2,3 ' 



,(111) _ 1 - h - b 2 + bi2 + 621 ~ frl21 

(1 + q + q 2 ){\ + q) 



r i,2,3 : - Cl , „ 1 „2\/1 I „1 ■ ^. iy - ) 



Therefore, by subtracting the full antisymmetrizer Yy 2 3 from /1, the first Young operator Y± 3 2 

(21) 
of a mixed type was found. Then, the second Young operator Yj 2 3 of a mixed type was computed 

by applying the reversion (conjugate) to Y x l 3 2 , that is, Y 1 \ 2 ' 3 := Y xa2 ~. NamelyQ 

> Y21 . 132 : =bexpand(f . 1-Y111 . 123) : Y21 . 123 : =bexpand(reversion(Y21 . 132) ) : 

Furthermore, Y^ 2 3 ' annihilates ly 3 2 when multiplied from both sides. This is a reflection of the 
fact that the left ideal (the representation space) in the Hecke algebra H$(3,q) generated by /1 
decomposes into a direct sum of one-dimensional left ideal and a two-dimensional left ideal. 

> simplify(Y.111.123+Y.21.132-fl) ; 



> simplify(Y.21.123+Y.3.123-reversion(f 1)) ; 



> CS(Y.111.123 &c Y.21.132),CS(Y.21.132 &c Y. 111. 123); 

0, 

> CS(Y.21.132 &c Y.3.123),CS(Y.3.123 &c Y. 21. 123); 

0, 

> CS(Y.111.123 &c Y. 111. 123 - Y. Ill . 123) ,CS(Y. 21 . 132 &c Y. 21. 132 - Y. 21. 132); 

0, 
In addition to the symmetrizer Y X23 and the antisymmetrizer Y± 2 3 displayed in rtlq) and (fly), 

( 2 l) f21) 

we have also two Young operators ly 2 3 and Y± 3 2 of mixed symmetry: 

> , Y21.123 , =bexpand(Y21.123); 

Y21 123 - qM + {Kiq3 + 2Kiq2 - q2 + 2Kiq + Ki) bl 

q+l + q 2 (q + l + q 2 )(l + q) 

_ (K 4 q 3 + 2 K 4 q 2 + q + 2 K 4 q + K 4 ) q b2 
(q+l + q 2 )(l + q) 
(K 4 q 3 + 2 K A q 2 + q + 2 K 4 q + K A ) b!2 

q 3 + 2 q 2 + 2 q + 1 
(q 5 K 4 + q 4 K 4 + q 3 + K 4 q 3 - q 2 + K 4 q 2 + K 4 q + K 4 - 1) b21 
+ q{q 3 + 2q 2 + 2q + l) 

{q 4 K 4 + K 4 q 3 + q 2 -K 4 q-K 4 + 1) b!21 
+ (q + l + q 2 )q(l + q) 

> 'Y21.132'=bexpand(Y21.132); 



7 Due to lengthy displays, we refer Reader to formulas (33) and (34) in 
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¥2112,2= qM (K 4 q 3 + 2K 4 q 2 + 2K 4 q-l + K 4 )bl 



q+l+q 2 (q + l + q 2 )(l + q) 

(g 4 K 4 + 2 K 4 q 3 + 2 K 4 q 2 + K 4 q + 1) b2 
+ (g + l + g 2 )(l + g) 

(K4, q 3 + 2 K 4 q 2 + 2 K 4 q - 1 + K 4 ) bl2 
+ q 3 + 2 q 2 + 2 q + 1 

(g 5 iT 4 + g 4 ^4 + ^4 q 3 + g 3 + ^4 q 2 + ^4 g + g + ^4 - 1) b21 

q{q 3 + 2q 2 + 2q + l) 
(q 4 K 4 + K 4 q 3 + q 2 -K 4 q-K 4 + 1) bl21 
(g + l + g 2 )g(l + g) 
In order to represent our Young operator Y 1 3 2 as a product of a row-symmetrizer i?(13) and 
a column-antisymmctrizcr C(12) = /i, we use /i denned above and compute i?(13) from the 
equation 

Y<£> = fl(13)/ x . (20) 

In order to solve the above equation in Maple for i?(13), we need to find an element Y in the Hcckc 
algebra which would not only satisfy equation (^) but also such that Y + Y~= 1. This is because we 
want the column antisymmetrizer C(12) to remain related to i?(12) through the reversion. Notice 
also that we are justified in defining C(12) as equal to /i (modulo a normalizing factor) because 
/i generalizes the antisymmetrizer C(12) from S2 to S3 : 

> bexpand(cl2/(l+q)) ; 

Id bl 

1+g 1+g 

> bexpand(f 1) ; 
Id Kihl+Kiqb2+Klbl2 (K 4 q 3 + q + K 4 -l)b21 (-K 4 + K 4 q 2 + l)bl21 



1 + g ! " g(l + g) g(l + g) 

That is, /1 is seen to contain C(12)/(l + q) if we replace K 4 with 1/(1 + q). Thus, we first express 
Y in the Hecke basis contained in the list bset below and then we make sure that Y + Y~= 1 : 

> bset: = [Id,bl,b2,bl2,b21,bl21] : 
Y:=bexpand(add(P[i]*bset [i] ,i=i. .6)) : 
Y:=bexpand(op(clisolve2(Y+reversion(Y)-Id,Y))) ; 

Y : = ( \Pi q - \P Q q + \ P 6 g 2 + \p 2 g - \ P2 + \ - \ P 3 ) Id + P 2 bl + P 3 b2 

+ (-P 5 -Pe + Pe q) blS + P 5 121 + P 6 b!21 
Next we require that the above found element Y satisfies equation (pOJ): 

> R.13:=bexpand(op(clisolve2(Y21.132-cmul(Y,f_l) ,Y))) ; 
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niq qld {-q 2 + q 2 P 3 + P 3 q-l + P 3 )bl (q 2 P 3 + P 3 q + P 3 - 1) b!2 
tild := ; 7TT r -T 3 06 H : ^- 

1 + g (g + l + g 2 )g q{q+l + q 2 ) 

(q 5 P 3 + g 4 P 3 + g 3 Ps ~ q 2 + g 2 P3 + P 3 g - 1 + fb) 6gJ 

(l + g)g 2 (g+l + g 2 ) 

(g 4 P 3 + g 3 P 3 ~ Pa, q + 1 - P 3 ) b!21 
q 2 (q 3 + 2q 2 + 2q+l) 
We verify that i?(13) is an idempotent: 

> CS(R13 &c R13-R13); 


It was pointed out in J7J that when the Clifford product i?(13)/i is computed, the free parameter 
P 3 disappears: 

> bexpand(cmul(R13,f 1)) ; 

qld (K 4 q 3 + 2 K A q 2 + 2K A q-l + K A )bl 



g + 1 + g 2 %1 

(g 4 Ki + 2 K 4 q 3 + 2 K 4 q 2 + K A q + 1) b2 
+ %1 

(K 4 g 3 + 2 K 4 q 2 + 2 K 4 q - 1 + if 4 ) &i^ 
+ %1 

(g 5 K 4 + q 4 K 4 + K 4 q 3 + q 3 + K 4 q 2 + K 4 q + q + K 4 - 1) b21 

g%l 
(g 4 X 4 + K 4 q 3 + q 2 - K 4 q - K 4 + 1) bl21 
(g + l + g 2 )g(l + g) 

%1 :=g 3 + 2g 2 + 2g+l 

(21) 
which gives z^ 3 2 computed above: 

> CS(7„-Y21.132); 



In order to construct the representation spaces from the four Young operators, one needs to find 

at least one Garnir G\ element in the Hecke algebra (7| |2Cj. All Garnir elements can be seen 

to act as row or column cycles in Young tableaux thereby generating non-standard tableaux which 

correspond to the basis vectors of the representation space. A Garnir element has the following 

defining properties: 

^SM 2 ! = 0- (21) 

G { ?i Y $l + 0. (22) 

The reason for requiring (|22| ) is that we want G\ 1 Y± 2 3 to be a second basis element in the left 

Hecke ideal generated by Y 1 ^ 2 3 . In order to solve (Elf) and (123) , we begin by assigning to X a 
general element of the Hecke algebra, that is, X is a linear combination of the Hecke basis elements 
with some undefined coefficients Ki, i — 1, . . . , 6 : 
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> X : =bexpand(K [1] *Id+K [2] *bl+K [3] *b2+K [4] *bl2+K [5] *b21+K [6] *bl21) ; 

X = K 1 Id + K 2 bl +K 3 b2+ K A bl2 + K 5 b21 + K 6 bl21 
We use clisolve2 to solve equation (|2l|). All three solutions returned by Maple are assigned to a 
list sol. 

> sol :=clisolve2(cmul(Y21. 123, X) , [seq(K[i] ,i=l. .6)] ) :nops(sol) ; 

3 

> for i from 1 to nops(sol) do X. i :=bexpand( subs (sol [i] ,X) ) od; 



XI := (K 2 q-K 6 q+ K 4 ) Id + K 2 bl +K 3 b2 + K 4 bl2 
(K 6 q 2 + K 4 q 2 -K 3 q-K 6 q-K 4 q + K 2 + K A ) 121 



K 6 b!21 



X2^K 1 Id + ^ 3 + ^ K ^ K6q \ +qKl + 1)bl + K 3 b2- M2 



q 2 (i + q) q(i + q) 

(q 5 K 6 -q 4 K 3 -q 3 K 3 -q 3 + K (> q 2 + q 2 K 1 + q 2 + qK 1 -q+l)b21 

T77" "1 — ^ r 6 ° Uli 

q d {l+q) 

X3 := K x Id 

(q-l + K 6 q + 2K 6 q 2 + 2K 6 q 3 + K 1 + 2qK 1 + q 4 K 6 + 2q 2 K 1 + q 3 K x ) bl 



q(q 3 + 2q 2 + 2q+l) 

K * w - , :, .,••"" ■ : + ((<? 6 k 6 + q 5 K 6 - q 5 K 3 -2 q * k 3 - g 4 + g 4 k 6 



(g-i)M* . ...,•■,■• ..-■,■■ ..-. 

q 3 + 2 q 2 + 2 q + 1 



+ 2 q 3 + q 3 K 1 +K 6 q 3 -2 q 3 K 3 + K 6 q 2 + 2 q 2 K x - q 2 K 3 - 2 q 2 + K 6 q + 2 qK x 

+ 2q + K 1 - 1)621 )/{q 2 (q 3 + 2q 2 + 2q+l)) + K 6 1121 
One way to find out if the three solutions returned by Maple are really linearly independent is to use 
a procedure f indbasis which from the given list of Clifford polynomials extracts all polynomials 
which are linearly independent. 

> nops (f indbasis ( [X . 1 , X . 2 , X . 3] ) ) ; 

3 
The three solutions X\, X 2 ,X 3 are therefore linearly independent. Another way to verify that fact 
would be to try to find three coefficients c±,c 2 ,c 3 , not all equal to zero, that would satisfy the 
following linear combination: 

c 1 X 1 +c 2 X 2 +c 3 X 3 = Q. (23) 

Thus, we can use again the procedure clisolve2 and try to solve equation (|23j) as follows: 

> clisolve2 (c [1] *X . 1+c [2] *X . 2+c [3] *X . 3 , [c [1] , c [2] , c [3] ] ) ; 

[{a = 0, c 2 = 0, c 3 = 0}] 
As expected, all coefficients are zero. Thus, Xi,X 2 ,X 3 are three linearly independent solutions of 
@, that is, Y^flXi = for i = 1, 2, 3. 

We proceed now to verify whether the solutions we have just obtained satisfy also equation (|2J 
In particular, we will try to see if there any non-zero values of the parameters K 2 , K 4 , K§, Kq in 

(21) 

Xi so that a i1y 2 3 would be 0. We will again use the procedure clisolve2. 

> varl : =select (type , indets (XI) , indexed) ; 
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varl := {K 2 , K±, K 5 , K 6 } 

> clisolve2(X.l &c Y. 21 . 123, varsl) ; 

[] 
As expected, Maple returns an empty solution set. 

3.3 Automorphism a q and the Garnir elements in the Hecke algebra 

In M, an automorphism a q was introduced in the Hecke algebra via the formula (48). a q re- 
places the reversion ~ and gives the inverse of the basis element bj< for any multi-index K. The 
automorphism a q is then extended to the whole Hecke algebra by the following definition: 

a q (b n ...b is ) = (-±y(b n ---b ls y 

= (^-) s (bl---b~ 1 ) (24) 

= a q (bi 3 )- ■ -a q (bi t ). 

In CLIFFORD, the a q automorphism has been programmed as a Maple procedure alpha2. For 
example, when s = 2, verification of (p3) can be done as follows :H 

> evalb (bexpand(alpria2(bl &c b2) )=bexpand(alpha2(b2) &c alpha2(bl) ) ) ; 
evalb (bexpand(alph.a2(b2 &c bl) )=bexpand(alpha2(bl) &c alpha2(b2))) ; 

true, true, 

The following is a verification that indeed ct q (bi) = bj~ = 1 + — , i = 1, 2.[] 

> bexpand(alpha2(bl) ) ,evalb(bexpand(cinv(bl) ) = bexpand(alpha2(bl) ) ) ; 
bexpand(alpha2(b2) ) ,evalb(bexpand(cinv(b2) ) = bexpand(alpha2(b2)) ) ; 

(q-l)Id bl 

H , true 



q q 

(q - 1) Id b2 



, true 



q q 

On the other hand, ( |2~4| ) implies that 0^(612) = (— ) 2 612 = b^ 

> bexpand(alpha2(bl2)) ; 

map (normal ,bexpand( (-1/q) ~2* (reversion (bl2)) ) ) ; 
evalb (bexpand (cinv (bl2) ) =bexpand (alpha2 (bl2) ) ) ; 



(1- 


-2q + q 2 )Id (q-l)bl (q - 1) b2 b21 


(1- 


q 2 q 2 q 2 q 2 
-2q + q 2 )Id (q-l)bl {q - 1) b2 121 

? 2 9 ' 2 

q z q A q z q A 



true 
and similarly for 621- For completeness we only show that 0^(6121) = ( — ) 3 &i2i = b~[ 21 : 

8 Procedure alpha2 is part of a package SUPPL. Rather than displaying and comparing long expressions, it is often 
convenient to use Maple's built-in Boolean procedure evalb which returns true or false when the two expressions 
are equal or not. 

9 In CLIFFORD, the symbolic inverse of any element can be found using a procedure cinv. 
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> bexpand(alpha2(bl21)) ; 

map (normal ,bexpand((-l/q) "3* (reversion (bl21) ) ) ) ; 
evalb(bexpand(cinv(bl21))=bexpand(alpha2(bl21))) ; 

(-2q 2 + q 3 + 2q-l)Id (l~2q + q 2 )bl {l-2q + q 2 )b2 (q - 1) b!2 

q 3 q 3 Q 3 I 3 

(q - 1) b21 bl21 

+ — - 1 

q 3 q 3 

{-2q 2 + q 3 + 2q- I) Id {l~2q + q 2 )bl {l-2q + q 2 )b2 (q - 1) b 12 
q 3 q 3 q 3 q 3 



(q - 1) b21 bl21 
q 3 q 3 



true 
■=1\\K\ C, - h-1 



Thus, we have verified that a g (6j<-) = ( — ) bn — b K for a multi-index K of length \K\. However, 
this property does not extend to non-homogeneous (non-versor like) elements of the Hecke algebra. 
Let hecke be a Maple variable representing an arbitrary element in -Hf(3, q) expanded in the Hecke 
basis {IMMM2M1M21} '■ 

> hecke :=h[l]*Id + h[2]*'bl> + h[3]* , b2 ) + h[4]*'bl2' + h[5]*'b21' + h[6]*'bl21'; 

hecke := h\ Id + h 2 bl + h 3 b2 + /i 4 bl2 + h 5 b21 + h e bl21 
Then the action of a q on hecke can be written as follows: 

> alpha2he eke : =bexpand ( alpha2 (he eke ) ) ; 

HI Id H2 bl H3 b2 H4 b!2 H5 b21 H6 bl21 



alpha2hecke :- 



q 3 q 3 q 3 q 3 q 3 q 3 



where the expressions Hi,i = 1, ... ,6, are q -polynomials with coefficients expressed in terms of 
hi,i = 1, . . . , 6.F] With a little experimentation it can be easily verified that, for example, a q does 
not give the inverse of the element 1 + b\ : 

> h[l] ,h[2] ,h[3] ,h[4] ,h[5] ,h[6] : =1 , 1 ,0,0,0,0 :bexpand(hecke) ; 

Id + 11 

> bexpand(CS(cinv(hecke) ) ) , bexpand(alpha2 (hecke) ) ; 

1 (q-2)Id 1 bl (2q-l)Id bl 

2 q-l + 2?-l' q + ~q~' 

We return now to the problem of finding the Garnir element G\ x that would satisfy equations 
( fH| ) and (|2^). Recall that equation ( pl| ) had three linearly independent solutions Xi,X 2 , X 3 . Notice, 

that X\ does not annihilate the Young operator Y X2 3 and a q {X\) does not annihilate the Young 

(21) 
operator Y 1 3 2 n0 matter what values are assigned to the parameters K 2 , K4, if 5, Kq : 

> clisolve2(X.l &c Y. 21 . 123, varsl) ,clisolve2(alpha2(X. 1) &c Y. 21 . 132, varsl) ; 

— ('21') (y]} 

Furthermore, the following shows that the elements X\Y X2 -s an d OL q (Xi)Y^ Z2 are linearly inde- 
pendent: 

10 In order to get the display for a q (hecke) shown above expressions Hi,i = 1, . . . , 6, have been defined in Maple 
as aliases. Sec Appendix 2 for their definitions. 
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> nops(findbasis([X.l &c Y. 21 . 123,alpha2(X. 1) &c Y. 21. 132])); 

2 
Since the six elements 

/v (3) v {2\) y _(2i) ry s v (2i) v (2i) v (m)-, 

\ x 1,2,3) J 1,2,3)^1 J 1,2,3) a Q\- A i) I l,3,2> J 1,3,2) J 1,2,3 J 

are linearly independent, 

> nops(findbasis([Y.3. 123, Y. 21. 123.X. 1 &c Y.21 . 123,alpha2(X. 1) &c Y. 21. 132, 
Y. 21. 132, Y. 111. 123])); 

6 

they may form a basis S for the Hecke algebra Hp(3,q). We define the Garnir element G\ x to be 

equal to Xi.f^ 

> alias (tl=K [6] *q~2+K [4] *q~2-K [5] *q-K [6] *q-K [4] *q+K [2] +K [4] ) : G21 . 1 . 1 : =bexpand(Xl) ; 

t1 b2 
G2111 := (K 9 q-K 6 q + KA Id + K 2 bl + + K 4 bl2 + K 5 b21 + K 6 bl21 

q 

> S:=[Y.3.123,Y.21.123,G21.1.1 &c Y.21. 123, 
alpha2(G21.1.1) &c Y. 21 . 132, Y. 21 . 132, Y. Ill . 123] : 

Procedure yexpand (See Appendix 2) is used to expand elements in the Hecke algebra Hf(3,q) in 
terms of the Young basis S. For example, 

> yexpand (Y3. 123) ; 

Y3.123 

> yexpand(&c(G21.1.1,Y21.123)); 

G21.lAkcY21.l23 

> yexpand(&c(alpha2(G21.1.1),Y21.132)); 

a2(G21. 1. 1) kcY21. 132 

Recall that the original basis in the Hecke algebra was: {1, bi, 62, &12, &21, ^121}- Each original basis 
element should be representable in terms of the Young basis S : 

> ' Id'=yexpand(Id) ; 

Id= Y3.123+ Y21. 123+ Y21. 132+ Ylll. 123 

> ' b l ,= y expand (bl) ; 

w2Y21. 123 qwl {G21 .1.1 kcY21 .123) 

01 — Yo.A-Zo -\- - -\- - - 

1 + q w3{l + q) 

q 2 (a2(G21 .1.1) kcY21 .132) qw5 Y21. 132 
w4 w6 

> ' b2'=y expand (b2 ) ; 



q Ylll. 123 



b2= Y3 123 - w7qY21A23 w8q{G21. 1.1 kcY21. 123) 



1 + q w3(l + q) 

q 3 (a2{G21. 1.1) kcY21. 132) w9 Y 21 .132 



w4 w6 



q Ylll. 123 



11 The above computations could have been performed with X2 or X3 instead of X\, and the results would have 
been similar. 
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where Wi,i = 1, . . . , 9, are polynomials in q parameterized in terms of Kq, K4, K5, Kq. They have 
been defined as Maple aliases and are shown in Appendix 2. 

4 Singular Value Decomposition 

Our next application of Clifford algebras will be to the Singular Value Decomposition (SVD) of a 
matrix (2^1 ■ There are many uses of SVD such as in image processing, description of the so called 
principal gains in a multivariable system fl25| , or in an automated data indexing known as Latent 
Semantic Indexing (or LSI). LSI presents a very interesting and useful technique in information 
retrieval models and it is based on the SVD ]l0|| . While in these practical cases computations 
are done numerically, it may be of interest to ask whether such decomposition of a matrix can 
be performed in the framework of Clifford algebras. That is, if any new insights, theoretical or 
otherwise, into such decomposition could be gained when stated in the Clifford algebra language. 
In this section we will present examples of such computations. 

We will explore a well-known fact that when p — q ^ 1 mod 4, Clifford algebra C£ p ^ q is a simple 
algebra of dimension 2™, n = p + q, isomorphic to a full matrix algebra Mat(2 fe ,K) of 2 fc x 2 k 
matricesQj with entries in K which is R, C, or H (see ||). Thus, any operation performed on a 
matrix A can be expressed as an operation on a corresponding to it element p in C£ p . q . The choice of 
the signature (p, q) depends on the size of A and the division ring K. Of course, for computational 
reasons one should find the smallest Clifford algebra C£ Piq such that the given matrix A can be 

embedded into Mat(2 fe , K) ~ C£ p>q . In the following we will use the same approach as in B where a 
technique for matrix exponentiation based on the isomorphism tp was presented. In particular, we 
will use a faithful spinor representation of C£ PA in a minimal left ideal S = C£ p . q f generated by 
a primitive idempotent /. Symbolic computations of such representations with CLIFFORD were 
shown in [0. 

Following p3], let A be an m x n real matrix of rank r. Then the SVD of A is defined a 
factorization of A into a product of three matrices U, E,!^ 1 where U and V are orthogonal 
matrices m x m and n x n respectively, and £ is a m x n matrix containing singular values of A 
on its "diagonal" . 

A = UY,V-\ U T U = I, V T V = I. (25) 

The matrices V = [wil^l • ■ ■ \ v n] an d U — [ui|w2| • ■ ■ \u m ] contain orthonormal bases for all four 
fundamental spaces of A. Namely, the first r columns v\, V2, ■ ■ ■ ,v r of V provide a basis for the 
row space TZ{A T ) while the remaining n — r columns of V provide a basis for the null space Af(A). 
Likewise, the first r columns ui,U2, • ■ ■ ,u r of U provide a basis for the column space C(A) while 
the remaining m — r columns of U provide a basis for the left-null space M{A T ). Vectors Vi are 
the normalized eigenvectors of A T A while vectors ut are the normalized eigenvectors of AA T . For 
i = 1, . . . , r, these vectors can be chosen to be related via the positive singular values <ii of A which 
are just the square roots of the eigenvalues of A T A (or of AA T .) Namely, 

Avi = aiUi, i=l,...,r. (26) 

It is a little tricky to make sure that the above relation is satisfied: this is because the choice of 
vectors u, is independent of the choice of vectors Vi. However, it is always possible to do so as we 



12 The value of k = q — r q — p , where r; is the Radon-Hurwitz number. The Radon-Hurwitz number is defined by 
a recursion as Ti+g = r; + 4 and these initial values: tq = 0, r± = 1, r2 = r$ = 2, 7*4 = 7*5 = r$ = ri = 3. 
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will see below (see also p8[). In order to complete the picture, the orthonormal set {v\, . . . ,v r } 
needs to be completed to a full orthonormal basis for R™ while {u\, . . . , u r } needs to be completed 
to a full orthonormal basis for R m . Since the additional vectors are being annihilated by A and A T 
respectively, that is, they are eigenvectors of A and A T (or of A T A and AA T ) that correspond to 
the eigenvalue 0, care has to be exercised when finding them. For example, while the eigenvectors 
of the symmetric matrix AA T are automatically orthogonal provided they correspond to different 
eigenvalues, eigenvectors of AA T that correspond to the eigenvalue don't need to be orthogonal: 
in this case the Gram-Schmidt orthogonalization process is used to complete the two sets. 

4.1 Singular Value Decomposition of a 2x2 matrix of rank 2 

In this section we present our first example of SVD applied to a 2x2 real matrix of rank 2. 
The purpose of this example is just to show step by step how finding the SVD of a matrix can be 
done in the Clifford algebra language. Reader is encouraged to perform these computations with 
CLIFFORD and an additional package ASVD which is described in Appendix 3. 
> A:=matrix(2,2, [2,3,1,2]) ;#defining A 

m:=rowdim(A) : #number of rows of A is m 

n: =coldim(A) : #number of columns of A is n 

2 3 



A:= 



1 2 



Since A £ Mat(2,R), we need to find (p, q) such that C£ p>q ~ Mat(2,R). As shown next, we have 
two choices for the signature: 

> all_sigs(2. . 2, real, simple) ; 

[[1, 1], [2, 0]] 
Thus, we can pick either Ct\ t \ or Cl%$. Our choice is C^2,o- We define a bilinear form B = diag(l, 1) 
and display information about C(.2,o- 

> dim: =2:B: =diag(l , 1) : eval(makealiases (dim) ) :data:=clidata() ; 

data := [real, 2, simple, — Id -\ — el, [Id, e2], [Id], [Id, e2]] 

The above output means that C£2,o is a simple algebra isomorphic with Mat(2, R); that the element 
h + hei is a primitive idempotent which we will call /; that the list [Id, e2] shown as the fourth 
entry displays generators of a minimal left-ideal C^2.o/ considered as vector space over R; that 
the division ring K = /C^2,o/ =< Id >r— R; and that the last list [Id, e2] gives generators of 
C(-2,of over K, and since K ~ R, it is the same as the fourth list.[j. In the following, we define 
a Grassmann basis in C£2.o, assign the primitive idempotent to /, and generate a spinor basis 
in C£ 2 ,of- 

> clibas :=cbasis(dim) ; #ordered basis in CI (2,0) 

clibas := [Id, el, e2, el2] 

> f :=data [4] :#a primitive idempotent in CI (2,0) 
SBgens :=data[5] :#generators for a real basis in S 
FBgens :=data[6] :#generators for K 

13 For more information see [d| and CLIFFORD'S help pages. 
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Here SBgens is a IK -basis for S = Ctzfif- Since for the signature (2,0) we have K ~ R, S ~ R 2 , 
and C^2,o — Mat(2, R), the output from the procedure spinorKbasis shown below has two basis 
elements and their generators modulo / : 

> Kbasis :=spinorKbasis (SBgens , f ,FBgens , 'left ' ) ; 

Kbasis := [[- Id + - el, - e2 - - el2], [Id, e2], left] 
Thus, the real spinor basis in S consists of the following two polynomials: 

> for i from 1 to nops (Kbasis [1] ) do f . i :=Kbasis [1] [i] od; 

fl :=-Id + -el, f2:=-e2-- el2. 
2 2 2 2 

Now, we compute matrices M\, M2, M3, M4 representing each of the four basis elements {1, ei, e2, e^} 
in (74,0-0 

> for i from 1 to nops(clibas) do M[i] :=subs (Id=l ,matKrepr (clibas [i] ) ) od: 

We will use a procedure phi which gives the isomorphism <p from Mat(2,R) to Ct^fi- This way we 
can find the image p in (74,0 of any real 2x2 real matrix A. Knowing the image (p(Mi) of each 
matrix Mi in terms of some Clifford polynomial in C^2,o, we can easily find the image p = ip(A) 
of A as follows:^ 

> p:=phi(A,M); #finding image of A in Cl(2,0) 

p:= 2Id + 2e2 + el2 

> pT:=phi(t(A) ,M) ; #finding image of t(A) in Cl(2,0) 

P T :=2Id + 2e2- el 2 
Next, we compute a symmetric matrix A T A (denoted in Maple as ATA), its characteristic polyno- 
mial, eigenvalues, and its orthonormal eigenvectors i>i,i>2- Vectors v\ and i>2 will become columns 
of an orthogonal matrix V needed for SVD of AQ: 

> ATA:=evalm(t(A) &* A); #finding matrix ATA 

5 8 
8 13 

> pTp:=phi(ATA,M) ; #finding image of ATA in Cl(2,0) 

pTp := 9 Id -4 el +8e2 
which should be the same as 

> , pTp , =cmul(pT,p) ; 

pTp = 9 Id -4 el +8e2 
The minimum polynomial of A T A (or pTp — (p(A T A)) is: 

> climinpoly (pTp) ; 

x 2 -18.T+1 

and, in this case, it is the same as the characteristic polynomial of A T A : 

> pol :=charpoly (ATA,x) ;#characteristic polynomial of ATA 

pol := x 2 - 18 x+ 1 

14 Since a similar computation was done in |4l B[, we won't display the matrices. 

15 From now on, as use Maple alias (t=transpose), that is, t(A) denotes the matrix transposition in Maple. 

16 We use procedure radsimplif y to simplify radicals in matrices and vectors. 
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ATA 



In order to find eigenvalues and eigenvectors of A T A, we will use Maple's procedure eigenvects 
modified by our own sorting via a new procedure assignL. The latter displays a list containing two 
lists: one has the eigenvalues while the second has the eigenvectors. PI In the following, we will assign 
the eigenvalues to Ai, A2 and the (un-normalized, but orthogonal) eigenvectors we assign to v\, v%. 
> P: =assignL(sort ( [eigenvects (ATA)] ,byeigenvals) ) ;N:=P [1] : 



P := [2, [9 + 4 v / 5, 9-4^5], [ 



1. 



I I 



y/l 



1. 



I I 



V5 



for i from 1 to N do lambda. i :=P [2] [i] ; v. i :=map (simplify .normalize (P [3] [i] ) ) od: 
We can now verify that vectors v\,v^ are eigenvectors of A T A with the eigenvalues Ai, A 2 . 

> for i from 1 to N do map(simplif y ,evalm(ATA &* v.i - lambda. i*v. i) ) od; 

[0, 0], [0, 0] 
Similar verification can be done in 6^2,0 since one can view the 1 -column eigenvectors v±, v% as 
one-column spinors in S. We simply convert the two vectors to spinors svi , SV2 which we express 
in the previously computed spinor basis /1, /b- 

> spinorbasis: = ["fl" , "f2"] : 

for i from 1 to N do sv. i :=convert (v. i , spinor ,spinorbasis) od; 



svl 



fl 



(l + V5)/8 



sv2 



fl 



(-l + V5)/8 



VTo + 275 \/io + 2V5 V10-2-V/5 V10-2-V/5 

Since V\,V2 are eigenvectors of ATA, spinors sv\, SV2 must be eigenspinors of pTp = ip(A T A). 

> for i from 1 to N do simplify ( (pTp - lambda. i) &c sv.i) od; 

0,0 
We are now in position to define the orthogonal matrix V — [i>i|i>2]- 

> V:=radsimplif y (augment (v. (1 . . N) )) ; #defining matrix V 

1 „ 1 



V :-- 



Vio+iTI 

l + Vb 



A/10-2V5 
-\/5 + l 



V10 + 2V5 \/lO-2V5 
Since later we will need images of V and V T under ip, we compute them now and store under 
the variables pV and pVt. The fact that V is orthogonal can be easily verified in the matrix 
language; in C^2,o it can be done as follows: 

> simplif y (cmul(pVt ,pV) ) ; 

Id 
Now we repeat the above steps and apply them to AA T . In the process, we will find its eigen- 
vectors Ui,U2. We must make sure that Av t = OiUi where c^ = VAj", i= 1, 2. This will require 
extra checking and possibly redefining of the u 's. 

> AAT:=evalm(A &* transpose (A) ) ; #computing AAT 

13 I 



AAT := 



8 5 



The image of AA T under ip in C^2,o w e denote as ppT. 

17 The first entry 2 in the output is just the number of eigenvectors. 
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> ppT:=phi(AAT,M) ; #finding image of AAT in Cl(2,0) 

ppT := 9 Id + 4 el + 8 e2 
In this case, the minimal polynomial of ppT and the characteristic polynomial of AA T are the 
same. 

> pol2: =charpoly( AAT, lambda) ; #finding the characteristic polynomial of AAT 

pol2 := A 2 -18A + 1 

> 'ppT'=climinpoly(ppT) ; 

ppT = x 2 -18x+l 
Since matrices A T A and AA T have the same characteristic polynomials, their eigenvalues will be 
the same. We define therefore the singular values o\ and 02 of A : 

> for i from 1 to N do sigma. i :=sqrt (lambda. i) od; 

crl:=V5 + 2, cr2:=\/5-2 
When we compute the eigenvectors u\, 112 of AA T , we will not necessarily have Avi = <7iUi, i = 1, 2. 
This is because the choice of u\, U2 is not consistent with the choice of Ui,f2- 

> P: =assignL (sort ( [eigenvects (AAT)] ,byeigenvals) ) : 

for i from 1 to N do lambda. i :=P [2] [i] ; u. i :=map(simplif y , normalize (P [3] [i] ) ) od: 

However, Av\ = o\U\ while Av2 — —a^ui : 

> radsimplif y(evalm(A &* vl-sigmal*ul) ) ; #this one checks out but 

[0, 0] 

> radsimplif y(evalm (A &* v2-sigma2*u2) ) ; #this one does not check out 
radsimplif y(evalm (A &* v2+sigma2*u2) ) ; #this one does check out 

14 -6 a/5 8-4V5 



V 10 - 2 a/5 \A0-2\/5 

[0,0] 
Notice that the set {^1,1*2} is orthonormal, but so is {ui,— U2}. Let's re-define U2 as —U2 and 
call it U22- For completeness we rename u\ as u\i : 

> ull:=evalm(ul) :u22 :=evalm(-u2) : 

In the Clifford algebra CH.2,0, we need to perform similar computations with vi,v%. The images 
i^(uii), (^(1*22) contained in the spinor ideal need to be found first. We call them su\ and SU2- 

> for i from 1 to N do su. i :=convert(u. i . i, spinor ,spinorbasis) od; 

-<-( 1 + ^ h2 , g su ,-(- 1 + ^^ ■> /' 



VlO + 275 VlO + 275 a/10-2a/5 ^10-2^ 

The verification of the condition (GO) in (7^2.0 looks as follows: 
> for i from 1 to N do simplify(p &c sv. i-sigma. i*su. i) od; 

0, 
Now we may define the orthogonal matrix U — [un|u22] and its image f(U) in Ci.2,0 which we 
callp[/:0 

18 For a later verification we will also need pUt = ip(U T ). Expressions %1 and %2 showing up in the Maple output 
for pU are just place holders for \f 10 + 2 \/5 and V10 — 2\/5 respectively as shown at the end of the display. 
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> U: =radsimplify(augment (ull ,u22) ) ; #defining matrix U 

1 + V5 -1 + V5 



U := 



- l 2- X 



VlO + 275 ^10-2^5 

pU : =phi (U , M) ; pUt : =phi (t (U) , M) : 



Hy %2V5 + l%l 



l%lV5)e2 + (^%2V5-i%l + ^%l>/5)ei2 



%1 := V10T2V5 



%2 := a/10 -2^ 
The fact that [/ is an orthogonal matrix can be easily now checked both in the matrix language 
and in the Clifford language: 
> radsimplif y (evalm(t (U) &* U));#U is an orthogonal matrix 



> simplify (pUt &c pU) ; 



1 
1 

Id 



Finally, we define matrix E using a procedure makediag. Recall p8j that E has the same dimensions 
as the original matrix A and that E T £, EE T are the diagonal forms of A T A and AA T respectively. 
In this example matrices E T E and EE T are the same since E is a square diagonal matrix. Normally 
these matrices are different although their nonzero "diagonal" entries are the same. 



A 1 A = VT, 1 T.V 1 , AA 1 = C/EE J U 1 , E 



0-1 

a 2 



E J E = EE J 



o\ 
erf 



(27) 



> Sigma :=makediag(m,n, [seq(sigma. i ,i=l . .N)] ) ; 

STS,SST:=evalm(t (Sigma) &* Sigma) ,evalm(Sigma &* t (Sigma)); 

V5 + 2 
V5-2 



STS, SST := 



£:= 
(V5 + 2) 2 







(V5 + 2) 2 

(a/5 -2) 



(\/5-2) 2 

> pSigma,pSTS,pSST:=phi(Sigma,M) ,phi(STS,M,FBgens) ,phi(SST,M) ; 

pSigma, pSTS, pSST := VEld + 2 el, 9 itf + 4\/5 el , 9 Id + 4^ el 
We should be able to verify in C^2,o the following two factorizations of AA T and ^4 T A 

A T A = V A E T EV" T 
AA T = [7££ T £/ T 

like this: 



(28) 
(29) 
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> evalb(pTp=simplify(pV &c pSTS &c pVt)), evalb(ppT=simplify (pU &c pSST &c pUt)); 

true, true 
We check the SVD of A, which is A = WEV T ,qj in the Clifford algebra language: 

> evalb(p=simplif y(pU &c pSigma &c pVt)), 

true, 
where 



' pU ' =pU ; 



P U = Id(-%lV5--%2 



±%2V5) + (^%lV5 + l%2 + ^%2V5)el 



+ (i%2V5+i%l-i%l>/5)e2 + (i%2>/5-i%l+i%l>/5)eiS 



%1 := V10 + 2V5 

%2 := ^10-2^5 
' pSigma ' =pSigma ; 



pSigma = v5 Id + 2 el 



> ' pVt ' =pVt ; 



1 



pVt = Id (--L %1 ^5 + I %2 - -L %2 V5) + (4 %l>/5 + 5 %2 - -J- %2 \/5) ei 

zU o 4(J zU o 4L) 



(_L %2 ,/5 + l%f + 1%1 V5) el> + (l%2,/5- 



f 



1 



%1 - — %1 V5) ei2 
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%1 := ^10-2^ 

%2 := yioTiTH 

4.2 Singular Value Decomposition of a 3x2 matrix of rank 2 

In this section we will show our second example of SVD applied to a non-square matrix. The matrix 
will need to be embedded first into an appropriate matrix algebra before its image can be found in 
a suitable Clifford algebra. 
> C :=matrix (3, 2, [3,0,0,-1,0,1]) ; 

m:=rowdim(C) : #number of rows of C is m 

n:=coldim(C) : #number of columns of C is n 



C := 



Since our matrix C is 3x2, we will embed it into Mat(4,K) ~ C£ PtQ where the signature (p,q) could 
be either (2,2) or (3,1). Since the symbolic spinor representation of Ctsi was already computed 
in la], we will work with the signature (3, 1). 
> dim:=4:B:=diag(l$3,-l$l) : eval(makealiases(dim) ) : 
19 The SVD of A is not unique: For example, A = (—U)Ti(—V T ) is another such factorization. 



" 3 


1 





-1 





1 
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Let's recall information about Clz,\ stored in CLIFFORD^]: 

> data: =clidata() ; 

data := [real, 4, simple, ' ' cmulQ\- Id + -el, - Id + - e34), [Id, e2, e3 , e23], [Id], 

[Id, e2, e3, e23}) 
We begin by defining a Grassmann basis in (% i, by assigning the fourth entry in the list to a 
primitive idempotent /, and by generating a spinor basis in S = Ct^f. 

> clibas :=cbasis(dim) : #ordered basis in Cl(3,l) 
f:=data[4]; #a primitive idempotent in Cl(3,l) 
SBgens :=data[5] ; #generators for a real basis in S 
N:=data[2]: #dimension of the spinor representation 

/ := cmulQ(± Id + l -el , l -Id + ^e34) 

SBgens := [Id, e2, e3, e23] 
Thus, the real spinor basis in S = Cl^.if consists of the following four polynomials: 

> for i from 1 to N do f . i :=SBgens [i] &c f od; 

fl := -Id + - e34 + - el + - el34 , f2:=-e2 + - e234 - - el2 - - el234 
4 44 4^' J 4 44 4 

1111 1111 

fS := - e3 + - el - - elS - - elA , fl := - e23 + - e24 + - el23 + - el24 

4 44 44 44 4 

We compute matrices Mi, . . . , Mig representing each basis element in C£$\. These matrices can be 

computed also with the help of a procedure matKrepr which, being less general than spinorKrepr 

used earlier, is also simpler to use. We won't display these matrices since they can be found in S. 

> for i from 1 to nops(clibas) do M[i] :=subs(Id=l , matKrepr (clibas [i] ) )od: 

Before we can use the procedure phi that realizes the isomorphism Mat(4,]R) ~ Cl?,,\, we need to 
embed matrix C into Mat (4, R). This is accomplished with a procedure embed. 

> A:=embed(C) ;ml :=rowdim(A) :nl :=coldim(A) : 

3 I) 

A:= 






-10 





1 









> p:=phi(A,M); #finding image of A (and C) in CI (3,1) 

pT:=phi(t(A) ,M) ; #finding image of AT (and CT) in Cl(3,l) 

p:= -Id + el - -e23 - -e24 +- e34 + - el23 + - el24 + el34 

P T := - Id + el + - e23 - - e2A + - e34 - - el23 + - el24 + el34 
2 4 424 4 

Next we compute A T A, AA T , their images in the Clifford algebra C^i under if, their eigen- 
values, and orthonormal eigenvectors Vi and Ui, % = 1, . . . ,4, which will become columns of two 

20 In the following display, ' cmulQ' '( » * Id+ ■= * el, ^ * Id+ ^ * e34) denotes an unevaluated product of two non- 
primitive idempotents i * Id+ ^ * el and ^ * Id+ ^ * e34 in Of3,i, and cmulQ is a name of a simplified version of 
the procedure cmul. While cmul gives the Clifford product in C£(B), cmulQ gives the Clifford product in C£(Q). 



2(. 



9 




9 








2 







1 


-1 





; 





-1 


1 

















orthogonal matrices V and U. 

> ATA,AAT:=evalm(t(A) &* A),evalm(A &* t(A)) 

ATA, AAT := 



> pTp,ppT:=phi(ATA,M) ,phi(AAT,M) ; #finding images of ATA and AAT in Cl(3,l) 

pTp, ppT := 

11 7 11 7 11 7 1 9 1 9 

— Id + - el + — e34 + - el34, — Id + - el + - e24 + - e34 - - el24 + - el34 
4 4 44 4 4 2424 

The characteristic polynomial of A T A and AA T is (x — 9)(x — 2)x 2 while the minimal polynomial 
of pTp and ppT is x{x — 2)(x — 9). 

> PI :=assignL(sort ( [eigenvects(ATA)] ,byeigenvals) ) ; 
P2:=assignL(sort ( [eigenvects(AAT)] ,byeigenvals) ) ; 
for i from 1 to N do lambda. i :=P1 [2] [i] ; 

v. i :=map( simplify, normalize (PI [3] [i] ) ) ; 
u. i : =map( simplify, normalize (P2 [3] [i] ) ) ; 
od: 

PI := [4, [9, 2, 0, 0], [[1, 0, 0, 0], [0, 1, 0, 0], [0, 0, 1, 0], [0, 0, 0, 1]]] 

P2 := [4, [9, 2, 0, 0], [[1, 0, 0, 0], [0, -1, 1, 0], [0, 1, 1, 0], [0, 0, 0, 1]]] 
The third list in PI and the third list in P2 contain vectors Vi and u% that will make up matrices 
V and U respectively. Clearly, since V is the 4x4 identity matrix, <p(V) = (p{V T ) = 1. In order 
to translate this matrix picture Clz ii we need to find spinors svi — <f(vi) and sui — tp(ui) in 
S = C£ 3 ,if. 

> spinorbasis^C'fl",' 'f 2" , "f 3" , "f4"] : 

> for i from 1 to N do sv. i :=convert (v. i, spinor ,spinorbasis) ; 

su. i :=convert (u. i, spinor ,spinorbasis) od; 

svl :=fl, sv2 :={2, sv3 :=f3, sv4 := f4 



-^V2f2 + ^V2f3, su3 



^V2f2 + ^V2f3, sv4 :=U 



sul :— fl , su2 

Spinors svi are eigenspinors of pTp, while sui are eigenspinors of ppT with the eigenvalues Ai = 9, 
A2 = 2, A3 = 0, A4 = 0, namely: 

> for i from 1 to N do simplify ( (pTp - lambda. i) fee sv.i); 

simplify ( (ppT - lambda. i) &c su.i) od; 

0, 0, 0, 0, 0, 0, 0, 

> V: =radsimplif y(augment (v. (1 . .N)) ) : #defining matrix V (identity matrix) 
pV,pVt:=phi(V,M),phi(t(V),M): #finding images of V and t(V) in Cl(3,l) 

We check that Avi — (JiUi where <7j = \f\i by verifying this fact in C^3,i.qj 

> for i from 1 to N do sigma. i :=sqrt (lambda. i) od; 

crl := 3, cr2 := y/2, ct3 := 0, cr4 := 



21 As before, we are following engineering practice |25| of considering square roots of zero eigenvalues of AA T and 
A T A also as the singular values of a matrix. While it is convenient to do so for computational reasons, it is not 
conceptually correct since by the definition singular values of a matrix are always positive. 
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> for i from 1 to N do simplify(p &c sv. i-sigma. i*su. i) od; 

0, 0, 0, 
Thus, we may now define an orthogonal matrix U — [ui\u2\v,3\u4,] and find its image f(U) and the 
image of its transpose in Cl-$.\ : 

> U:=radsimplify (augment (seq(u. i,i=l . .N)) ) ; ttdefining matrix U 

10 

~y/2 \V2 

ly/2 \V2 



U := 















1 



> pU:=phi(U,M);pUt:=phi(t(U),M); 



pU :=-Id+-el -jV2e24-jV2e34+jV2 e!24 +jV? el34 



P Ut 



1 



1 



„Id + ^el - -V2e24- -V2e34 + -V2 el24 + - \/2 el34 

2 2 4 4 4 4 



The fact that U is orthogonal is reflected in the following: 

> pUt &c pU; 

Id 
Finally, it is just enough to find matrix S and verify SVD for tp(A) in the Clifford algebra Cl^,i- 

> Sigma: =makediag(ml ,nl , [seq(sigma. i , i=l . . N)] ) ;pSigma:=phi (Sigma, M) ; 

" 3 " 

y/2 







£ : = 



pSigma := Id (^ + \ y/% + (- \ V2 + |) el + (| + ± >/2) e^ + (-|v^+|) e^ 

> evalb(p=pU &c pSigma &c pVt) ; #SVD of phi (A) 

true 
Since the original matrix C was 3 by 2 and not 4 by 4, in order to find SVD of C we need 
to project out certain columns and rows out of the matrices U, £, and V. This can also be done 
internally in the Clifford algebra. The original matrix C has therefore this factorization: 

> U2,Sigma2, V2t :=submatrix(U, 1 . . m, 1 . .m) , submatrix (Sigma, 1 . .m, 1 . .n) , 
submatr ix (t ( V) , 1 . . n , 1 . . n) ; 



U2, £2, V2t := 
> evalm(C)=radsimplify(evalm(U2 &* Sigma2 &* V2t)); 
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4.3 Additional comments 

In this section we have shown that it is possible to translate the matrix algebra picture of the 
Singular Value Decomposition of a matrix A into the Clifford algebra language. Although we 
have not abandoned entirely the linear algebra formalism in our examples, e.g., we have computed 
the eigenvalues and the eigenvectors of A T A and AA T , and only then we have found images of 
the eigenvectors in the spinor space C£ p , q f, these computations including solving the eigenvalue 
problems can be done entirely in Ci Viq without using matrices. For example, if we consider element 
pTp — ip(A T A) in Cls.i from the last example, we can find its eigenvalues from its minimal 
polynomial while its eigenvectors can be found by solving the eigenvalue equation. For example, we 
know that the minimal polynomial of p Tp is 

> pol : =f actor (climinpoly (pTp) ) ; 

pol :— x (x — 2) (x — 9) 
hence the eigenvalues of pTp are 9,2,0. Of course, the minimal polynomial doesn't give us their 
geometric multiplicities. However, we can find them by solving the eigenvalue equation directly in 
C^3.i for each of the eigenvalues. Let's assign these three known eigenvalues to Ai, A2, A3 : 

> Iambdal,lambda2,lambda3:=9,2,0; 

Al, A2, A3:=9, 2, 
Let ip € C^3 ; i/ be an arbitrary spinor expressed in the spinor basis {/1, / 2 , /3, fi\ computed earlier, 
and let c±, . . . , C4 be its undefined (real) coefficients. 

> psi : =c [1] * 'i 1 > +c [2] * 'f 2' +c [3] * 'f 3' +c [4] * ' f 4' ; 

iP:= Cl fl+c 2 f2 + c 3 f3 + c 4 f4 
We will now use the procedure clisolve2 from the SUPPL package (see Section 2) to solve the 
eigenvalue equation 

pTpi;^X^. (30) 

First we solve it when A = Ai =9. 

> soll:=clisolve2(pTp &c psi - lambdal*psi , [c [1] ,c [2] ,c [3] ,c [4] ] ) ; 

soil := [{c 2 = 0, c 3 = 0, c 4 = 0, ci = ci}] 

> psil :=c [1] *'fl ' ; #eigenspinor of pTp with eigenvalue lambdal=9 

rl>l:=cifl 

Thus, the first solution to ( p0|) for Ai = 9 is a one-parameter solution that belongs to a one- 
dimensional subspace spanned by f\. That is, the geometric multiplicity of the eigenvalue Ai = 9 
is 1. Similarly for A = A2 = 2 : 

> sol2:=clisolve2(pTp &c psi - lambda2*psi , [c [1] ,c [2] ,c [3] ,c [4] ] ) ; 

sol2 := [{c 3 = 0, c 4 = 0, ci = 0, c 2 = c 2 }] 

> psi2: =c [2] * 'f 2' ; #eigenspinor of pTp with eigenvalue lambda2=2 

ip2:=c 2 f2 
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The second solution to (BOJ) for A2 = 2 is also a one-parameter solution that belongs to a one- 
dimensional subspace spanned by / 2 . The geometric multiplicity of the eigenvalue A 2 — 2 is also 1. 

> sol34:=clisolve2(pTp &c psi - lambda3*psi, [c [1] ,c [2] , c [3] , c [4]] ) ; 

sol34 := [{c 3 = c 3 , c 4 = c 4 , c 2 = 0, a = 0}] 
The third solution to ( pOJ ) for A = A3 = is parameterized by two parameters c 3 and c 4 and 
belongs therefore to a two-dimensional subspace of S spanned by {/3,/i}. This implies that the 
geometric multiplicity of A = A3 = is two. 

> psi3:=c [3] *'f 3' ; #eigenspinor of pTp with eigenvalue lambda3=0 

ip3:=c 3 f3 

> psi4:=c [4] *'f 4' ; #eigenspinor of pTp with eigenvalue lambda3=0 

^4:=c 4 # 
So, up to a normalizing scalar, the eigenspinors ipi , ip2 > ^3 > V^ gi ye previously computed normalized 
spinors sv\, sv 2 , sv 3 , st> 4 . Likewise for the image ppT of AA T : 

> soll:=clisolve2(ppT &c psi - lambdal*psi, [c [1] ,c [2] ,c [3] ,c [4] ] ) ; 

soil := [{c 2 = 0, c 3 = 0, c 4 = 0, c\ = ci}] 

> phil :=c [1] *'f 1 ' ; #eigenspinor of ppT with eigenvalue lambdal=9 

<t>l~dfl 

The first one-parameter eigenspace of ppT corresponding to the eigenvalue Ai = 9 is one-dimensional 
and it is spanned by /1. 

> sol2:=clisolve2(ppT &c psi - lambda2*psi, [c [1] ,c [2] ,c [3] ,c [4] ] ) ; 

sol2 := [{c 3 = c 3 , c 4 = 0, ci = 0, c 2 = -c 3 }] 

> phi2 :=-c [3] *'f 2'+c [3] * 'f 3' ; #eigenspinor of ppT with eigenvalue lambda2=2 

02:=-c 3 f2+c 3 f3 
The second one-parameter eigenspace of ppT corresponding to the eigenvalue A 2 = is also one- 
dimensional and it is spanned by {/ 2 , f 3 }. 

> sol34:=clisolve2(ppT &c psi - lambda3*psi, [c [1] ,c [2] , c [3] , c [4]] ) ; 

sol34 := [{c 3 = c 3 , €4 = c 4 , ci = 0, c 2 = c 3 }] 
The third solution to ( |30| ) is parameterized by two parameters c 3 and c 4 . Thus, we get a two 
dimensional vectors space spanned by {/3, / 4 }. 

> phi3:=c [3] *'f 2'+c [3] *'f 3' ; #eigenspinor of ppT with eigenvalue lambda3=0 

03:=c 3 f2 + c 3 f3 

> phi4:=c [4] *'f 4' ; #eigenspinor of ppT with eigenvalue lambda3=0 

04:=c 4 /^ 
So, again up to a normalizing scalar, <pi, 2 , 4> 3l </> 4 give the eigenspinors su\, su 2 , su 3 , su 4 computed 
earlier. 

5 Clifford algebras in robotics 

Clifford algebras Cl(V, Q) on a quadratic space (V, Q) endowed with a degenerate quadratic form 
Q and associated groups Spin, Pin, Clifford, etc., were studied in |3[ g] and [g, |lj, Q. In contrast 
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to the Clifford algebras of a non-degenerate quadratic form, these algebras possess a non-trivial 
two-sided nilpotent ideal called Jacobson radical. The Jacobson radical J is generated by the null- 
vectors in V which are orthogonal to the entire space V (that is, J is generated by the orthogonal 
complement V 1 - of V) . It is known (l^, fii|] that J contains every nilpotent left and right ideal in 
C£(V, Q). From the point of view of the spinorial representation theory of Clifford algebras used in 
Section EL an important difference is that C£(V, Q) does not possess faithful matrix representation 
when Q is degenerate. 

Let V — V'-LV 1 - where V is endowed with a non-degenerate part Q' of Q of signature (p, q). 
Let dim(V^) = d, hence p + q + d = dim (V). Let's denote C£(Q) as C£d, P , q - Then we have a direct 
sum decomposition C£d, P , q = G£ Ptq ® J into C£ Ptq -modules. It was shown in M that when d = 1 
this decomposition is responsible for a semi-direct product structure of the group of units C£\ pq 
of C£i tPt q and of all of its subgroups such as the Clifford group r(l,p, fc) and the special Clifford 
groups T ± (l,p,k) = T(l,p,k)nC£ lpq , where C£f pq (resp. Ctf pq ) denotes the even (resp. odd) 
part of Ctx,p,q- The Clifford group was defined as T(l,p,k) = {g G C£* lpq \gvg~ 1 G V,v G V}, 
that is, without a twist. Q Let N : F(l,p, k) — > C£* k be defined as N{g) = gg. Then we define 
the reduced Clifford groups as Fg(l,p, k) = ker N n F ± (l,p, k). The Pin(l,p, q) and Spin(l,p, k) 
groups are then: 

Pin(l,p, q) = {g G r(l,p, q) \ N(g) = ±1}, Spin(l,p, q) = {g G F+(l,p, q) \ N(g) = ±1}. (31) 

In preparation for our computations below, from now on we assume that p + q is an odd positive 
integer. Let G = {1 + vei\v G V ,e\ = 0} be a subgroup of C£\ pq . Then it was proven in [Q that 

Pin(l,p, 9 )=GxF±(p,fc), Spin(l,p,(7)=GxSpin(p,fc). (32) 

In the above, symbol xi denotes a semi-direct product with the group on the right acting on the 
group on the left. For example, it will be of interest to us to note that the homogeneous Galilei group 
of rigid motions Gq — R 3 x SO(3) in R 3 is isomorphic to SO + (1,0, 3) and it is doubly covered 
by Spin + (1, 0, 3), the identity component of Spin(l,0,3). For a similar result to (B2h when one 
considers the twisted Clifford group Fq,(1,p, k) and a twisted map N a : r Q (l,p, k) — > C£\ k defined 
as N a (g) — gg where ~ denotes the conjugation in Cl\ lViq , see |1^, |2^] . 

In the following two sections we will use approach and notation from Selig pTJ where the author 
denotes the degenerate Clifford algebra C£d,p,q as C(p 7 q, d). Furthermore Selig uses twisted groups 
and defines the Pin and Spin groups as follows: 

Pin(n) = {g G C(0, n, 0) : gg* = 1 and a(g)xg* G V for all x G V}, (33) 

Spin(n) = {g G C+(0, n, 0) : gg* = 1 and gxg* G V for all x G V}, (34) 

with * denoting the conjugation ~ in the Clifford algebra C{p,q,d). It is implicit in the definitions 
above that the actions of Pin and Spin on V are x i— > a(g)xg* and x i— > gxg* respectively. 

5.1 Group Pin(3) 

In this section we will perform some computations with Pin(3). In particular, we will find all 
possible forms of the elements in Pin(3) and verify some facts about that group. We begin by 

22 See Crumeyrolle |L4] for a definition of the Clifford group with the twist given by a : in Crumeyrolle's notation 
a denotes the principal automorphism or the grade involution in C£(Q). Then the twisted Clifford group is defined 
as r a (l, P ,k) = {g£Cei pq \a(g)vg- 1 &V,veV}. 
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assigning a diagonal matrix to the bilinear form B. Grassmann basis for C£o,3 will be stored in the 
variable clibas. Following Selig we re-name Clifford conjugation as a procedure star and define 
a Euclidean norm on V = M 3 as a procedure Enorm. We will also define some additional Maple 
procedures that will be useful below. 

> B: =diag(-l$3) ;eval(makealiases(3) ) : clibas :=cbasis(3) ; 

" -1 

B := 0-10 

0-1 

clibas := [Id, el, e2, e3, el2, el3, e23, el23] 

> star :=proc(x) conjugation(x) end: #star (conjugation) operation in Cl(0,3) 
Enorm :=v->simplif y(scalarpart (v &c star(v))): #(pseudo)Euclidean norm in V 
alpha :=proc(x) gradeinv(x) end: #alpha (grade involution) operation in CI (0,3) 
scalarprod:=(x,y)->scalarpart (l/2*(x &c star(y) + star(y) &c x)): 

> Pin_action:=(x,g)->clicollect (simplif y(alpha(g) &c x &c star (g) ) ) ;#action of Pin(3) 

Pin_action := (x, g) — > clicollect (simplify ((a (5) '&c' x) '&c' star(g))) 

> Spin_action:=(x,g)->clicollect (simplif y(g &c x &c star(g) ) ) ;#action of Spin(3) 

Spin_action := (x, g) — > simplify((g '&c' x) '&c' star(g)) 
Let v, Vi,V2 be three arbitrary vectors in R 3 with some undetermined coefficients expressed in a 
pseudo-orthonormal basis {ei, 62,63} : 

> V:=cl*el+c2*e2+c3*e3:vl:=cll*el+cl2*e2+cl3*e3:v2:=c21*el+c22*e2+c23*e3: 
Then the Euclidean norm in R 3 is: 

> Enorm (v) ; 

cl 2 + c2 2 + c3 2 
The action of Pin on C^o,3 is realized as the procedure Pin_action defined above. Let's verify 
Selig's claim (|27]], page 153) that when x, g are both in V, then gxg* automatically belongs 
to V : 

> Pin_action(v, vl) ; 

(ell 2 c3 - cl3 2 c3 - 2 cl3 cl2 c2 + c!2 2 c3 - 2 cl3 ell cl) c3 

+ (-ell 2 cl - 2 ell c!3 c3 + c!2 2 cl + cl3 2 cl - 2 ell cl2 c2) el 

+ (ell 2 c2 - 2 cl2 ell cl -2 cl2 cl3 c3 + cl3 2 c2 - cl2 2 c2) e2 
As we can see from the above, the output of Pin_action(v,vl) belongs to V. In order to check 
that indeed the action of Pin in V preserves the scalar product, we will first find all possible forms 
of g € Pin(3). Recall that according to ( [33| ) any element g G Pin(n) must satisfy two conditions: 
(1) gg* = 1 and (2) a(g)vg* e V for any v € V. Suppose that g is an arbitrary element in C£o,3,o 
expressed in CLIFFORD in terms of the Grassmann basis {1, ei, e2, ei2}[j 

> g:=add(x.i * clibas [i] , i=l . .nops (clibas) ) ; #a general element in Cl(0,3,0) 

g := xl Id + x2el + x3 e2 + x4 e3 + x5 el 2 + x6 el3 + x7 e23 + x8 el23 
We will now attempt to find conditions that the coefficients Xi,i = 1, ... ,8, must satisfy so that 
gg* = 1. We will again use the command clisolve2. In order to shorten its outputs, additional 
aliases Kj,j = 1, ... ,5, need to be defined (see Appendix 2). 

23 R,ecall that ei2 was defined above as an alias of ei A e2 with the command makealiases. 
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The first condition (1) gives: 

> sol:=clisolve2(cmul(g,star(g) )-Id, [x. (1 . .8)] ) ; 

t rr 4 — x7k5 — x5x4+x6x3 , ~ ^ ,- « 

so( := [{xi = , x8 = x8, x4 = x4, x7 = x7, xo = xo, xo = xo, xz = reo, 

X<5 

k4 —x5x4+x6x3, 

x3 = x3\, \x8 = 0, xl = — , x4 = x4, x7 = x7, x6 = x6, x5 = xo, x3 = x3, x2 = \, 

i,i x7 x7 

r^ -I rp J rp K 

{x7 = 0, x8 = 0, xl = —-, x3 — — — , x2 = x2, x4 = x4, x6 = x6, x5 = x5\, 
xo xo 

{x4 = 0, x7 = 0, x6 = 0, x8 = 0, x2 = x2, x5 = x5, x3 = x3, xl = k2}, 

{x7 = 0, x6 = 0, x5 = 0, x8 = 0, x2 = x2, x4 = x4, x3 = x3, xl = k1}] 

Thus, there are five different possible solutions, three of which requiring respectively that xq, xj and 
^8 be non-zero. Let's substitute these solutions into g. 

> for i from 1 to nops(sol) do g. i :=subs(sol [i] ,g) od; 

gl := _ (- x7K Z- x5x 4+ x6x3 )Id +K5el+ x3 e2 + x4e3 + x5 el2 + x6 el3 + x7e23 

x8 
+ x8 e!23 

k4H (-x5 x4 + x6 xS) el 

g2 := + ^ — + x3 e2 + x4 e3 + x5 el2 + x6 el 8 + x7 e23 

x7 x7 

nild „ . x4 x5 e2 , „ _ . n „ jn 

g3 := — — + x2 el + -?—— h x4 e3 + x5 e!2 + x6 e!3 

xo xo 

g4 := k2 Id + x2 el + x3 e2 + x5 el2 

g5 := k1 Id + x2 el + x3 e2 + x4 e3 

The above are five different types of g in C£o,3,o satisfying gg* = 1. 

> for i from 1 to nops(sol) do simplif y (cmul(g. i ,star (g. i) ) ) od; 

Id, Id, Id, Id, Id 
We need to make sure now that each gi displayed above satisfies also the second condition (2), 
namely, a(g)vg* is in V for any v 6 V. We begin with the simplest element g§. By computing 
the Pin group action on v and requiring that the result be a 1 -vector, we get for g$ : 

> Pin_action(v,g5) ; 

(2 «1 c3 x4 + 2 k1 cl x2 + 2 nl c2 x3) Id + (-2 x3 x4 c2 ~ 2 x4 x2cl -2x4 2 c3 + c3) e3 

+ (-2 x4 x2c3 - 2 x3 x2 c2 - 2 x2 2 cl + cl) el 

+ (-2 x3 2 c2 -2x4x3 c3 - 2 x2 x3 cl + c2) e2 
It should be clear from the above that since the coefficient of Id must be zero for any c\,ci,c-s, 
either X2 = £3 = £4 = or k\ — 0. Let e — ±1-Fj Thus, the former gives g = ±1, 

> g.5.1:=subs({x2=0,x3=0,x4=0},g5) ; 

g51 := eps Id 
while the latter gives 

24 In Maple one way to make e = ±1 is to define alias (eps=Root0f (_Z"2-1)) :. In the following, Maple outputs 
will contain the alias eps. 
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> g.5.2:=subs({kappal=0,x4=lambdal},g5) ; 

g52 := Al e3 + x2 el + x3 e2 
where Ai = ±y/l — x 2 — x\^\ We will collect all Pin group elements in a set Pin_group. 

> Pin_group:={g.5.1,g.5.2}; 

Pin.group := {eps Id, Al e3 + x2 el + x3 e2} 
Similarly, we consider 34. We assign the identity coefficient of the Pin action a(g)vg* to a 
variable eq and find a solution to the resulting two equations that will be parameterized by ci , C2 : 

> a: =Pin_action(v,g4) ; 

a := (2 x3 k2 c2 - 2 x5 x2c2 + 2 x5 x3 cl + 2 x2 k2 cl) Id + c3 e3 
+ (-2 x2 2 cl - 2 x5 k2 c2 - 2 x3 x2 c2 + cl - 2 x5 2 cl ) el 
+ (-2 x5 2 c2 -2x2x3 cl + 2 k2 x5 cl - 2 x3 2 c2 + c2) e2 

> eq: =collect (coef f (a, Id) ,{cl,c2}) :eql :=coef f (eq,cl) : eq2 :=coef f (eq,c2) : 
sol : = [solve ({eql , eq2},{x2,x5,x3})] ; 

sol := [{x5 = 0, x2 = \2, x3 = x3}, {x3 = 0, x2 = 0, x5 = x5}} 
In the above, A2 — ±\/l — x 2 . Likewise, we set A3 = zfc^/1 — x^.p°|The two new elements we assign 
to 541,5142 and add to Pin_group. 

> for i from 1 to nops(sol) do g.4. i :=simplif y (subs(sol [i] ,g4) ) od: 
Pin_group : =Pin_group union {g41,g42}; 

Pin_group := {eps Id, Al e3 + x2 el + x3 e2, A2 el + x3 e2, A3 Id + x5 el2} 
In order to continue with g 3 displayed above, we must make the assumption x e ^ known to 
Maple. Then, the action of 53 on a vector can be computed.^] 

> assume (x6>0,x6<0) ; a: =Pin_action(v,g3) : 

As in the previous two case, the quantity a is spanned by {Id,e\,e2,e^\. We will isolate the 
coefficient of the identity element in a, assign it to a variable eq, and then determine for which 
values of X2,Xi,x$,XQ it will be automatically zero for every choice of ci, 02,03. This will require 
solving a set of three equations {eqi,eq2,eqs} for X2, Xi,x§,XQ. Maple reminds us that xq ^ by 
displaying it as x6~ : 

> cliterms(a) ; 

{Id, e3, el, e2} 

> eq:=collect(coeff (a, Id) ,{c. (1. .3)}) ; 

„ (x6~ 3 x2- \4x4x6~)c3 „ (~\4x2 x6~ - x5 2 x4 x6~ - x4 x6~ 3 ) cl 

eq := —2 ^ 2 7. 

x6~ 2 x6~ 2 

(-A4 x5 x4 + x2 x5 x6~ 2 ) c2 

xlr 2 

> for i from 1 to 3 do eq. i :=coef f (eq, c . i) od: 
sol :=solve({eq. (1 . .3)} ,{x6,x4,x5,x2}) ; 

sol := {x2 = 0, x6~ = x6~ , x5 — x5, x4 = 0} 

25 In Maple we define alias (lambdal=RootOf (_Z"2+x2"2+x3"2-l)) : . 

26 Both are defined as aliases: alias (lambda2=Root0f(_Z"2-l+x3"2)) : alias (lambda3=Root0f (_Z~2-l+x5"2)) :. 

27 We won't display it due to its length. 



34 



This time we only have one solution which we then substitute into g$. 

> g31 : =subs(sol ,g3) ; 

g31 := —— + x5 el2 + x6~ el3 
xo 
where A5 is another alias (see Appendix 2). 

> Pin_group: =Pin_group union {g31}; 

X5Id 

Pin-group := {A3 Id + x5 e!2, —— + x5 el2 + x6~ el3, Al e3 + x2 el + x3 e2, 

X2el +x3 e2, eps Id} 
By continuing in the similar fashion with the elements gi and g\ , one can find all general types 
of Pin(3). Finally, all general elements of Pin(3) can be displayed: 

> ' Pin_group'=Pin_group; 

Pin.group := {A3 Id + x5 el2, — — + x5 el2 + x6~ el3, 

+ x5 el2 + x6 el3 + xT e23, A9 el + x3 e2 + x4 e3 + x8~ el23, 

x7 

Al e3 + x2 el + x3 e2, \2 el + x3 e2, eps Id} 
where A7 and Ag are displayed in the Appendix 2. It is a simple matter now to verify that all ele- 
ments of Pin displayed in Pin-group satisfy both conditions (1) and (2) from the definition (|33|). 

> for g in Pin_group do evalb(simplif y(cmul(g, star (g) )=Id) ) od;#Condition (1) 

true, true, true, true, true, true 

> for g in Pin_group do 

evalb(Pin_action(v,g)=vectorpart (Pin_action(v,g) , 1) ) od;#Condition (2) 

true, true, true, true, true, true 
We are now in position to verify Selig's claim J27J, page 153, that the scalar product on V = M 3 
defined in C£d, p , q as 

vi • v 2 = -(viVj + v 2 Vi), for any v i; v 2 G M 3 , 

is preserved under the action of the Pin group. Let Vi,V2 be two arbitrary 1 -vectors defined 
earlier. Procedure scalarprod that gives the scalar product may be defined as follows: 

> scalarprod :=(x,y)->scalarpart (l/2*(x &c star(y) + star(y) &c x)): 
for g in Pin_group do 

simplify (scalarprod (Pin_act ion (vl ,g) ,Pin_action(v2,g) ) -scalarprod (vl, v2) ) 
od; 

0, 0, 0, 0, 0, 
Thus, Pin(3) preserves the scalar product in R 3 and, therefore, we have a homomorphism from 
Pin(3) to 0(3) which is known to be a double-covering map. In the process, we have found all 
types of elements in Pin(3). 



:-!. r ) 



5.2 Group Spin(3) 

In this section we will perform a few computations with Spin(3). Once we have found general 
elements in Pin(3), it is much easier to find elements in Spin(3). Recall from ( pMJ ) that 

Spin(3) = {g£ C£^ 30 : gg* = 1 and gxg* E M 3 for all x e R 3 }. 

Let's find gSpin, a general element in Spin(3). Since Spin(3) C C£q 30 , we will begin with 
decomposing gSpin over the even basis elements. 

> clibaseven: =cbasis (3, 'even' ) ; 

clibaseven := [Id, el2, el3, e23] 

> gSpin:=c0*Id+c3*el2+c2*el3+cl*e23; 

gSpin := cO Id + c3 el2 + c2 el 3 + cl e23 
Notice that under the Spin group action defined as a procedure Spin_action 

> Spin_action: =(x,g)-> simplify (g &c x &c star(g)); 

Spin_action :— (x, g) — > simplify ((g '&c' x) '&c' star(g)) 
vectors are automatically mapped into vectors: 

> Spin_action(v,gSpin)-vectorpart (Spin_action(v,gSpin) , 1) ; 


We just need to make sure that gg* = 1 for each g S Spin(3). To simplify Maple output, we define 
K = y/l — cf — c| — c§ as a Maple alias. 

> alias (kappa=sqrt(-cl"2-c2"2-c3"2+D) : 

sol :=clisolve2(cmul(gSpin,star (gSpin) )-Id, [c . (0 . . 3)] ) ; 

sol := [{cl = cl, c2 = c2, c3 = c3, cO = k}, {cl = cl, c2 = c2, c3 = c3, cO = -k}] 

> gSpin:=eps*kappa*Id+c3*el2+c2*el3+cl*e23; 

gSpin := eps k Id + c3 el2 + c2 el3 + cl e23 
Thus, the most general element in Spin(3) is just g = enl + c^e^ + c^eis + cie23 where e = ±1. 
Notice, that the defining properties of g are easily checked: 

> simplif y (cmul(gSpin, star (gSpin) ) ) ; 

Id 

> evalb(Spin_action(v,gSpin)=vectorpart (Spin_action(v,gSpin) , 1) ) ; 

true 
In fact, element g G Spin(3) could be identified with a unit quaternion spanned over the basis 
{I,ei2,ei3,e23}. Then, the * conjugation becomes the quaternionic conjugation. It can be easily 
checked by hand or with CLIFFORD that the basis (bi)vectors anticommute and square to — 1. 

> quatbasis : = [el2,el3,e23] ; 

quatbasis := [el 2, el3, e23] 

> M:=matrix(3,3, (i,j)->cmul (quatbasis [i] .quatbasis [j] ) ) ; 

-Id c23 -el3 

M := -c23 -Id el2 
el 3 -el2 -Id 
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We have unit quaternions on a unit sphere in R 4 isomorphic to Spin(3) while the even part of C^o,3 
is isomorphic with the quaternionic division ring H. Spin(3) acts on R 3 through the rotations. 
In Appendix 3 one can find a procedure rot which takes as its first argument a vector v in R 3 
while as its second argument it takes a quaternion. For example, a counter-clockwise rotation in the 
{ei,e2J is accomplished with a help of a unit quaternion cos(|) + sin(|) e 12 : 

> rot(el,cos(theta/2)+sin(theta/2)*el2) : 
rot (e2 , cos (theta/2) +sin(theta/2) *el2) ; 
rot (e3 , cos (theta/2) +sin(theta/2) *el2) ; 

cos(fl) el + e2 sin(0), -el sin(0) + cos(9) e2, e3 
Let's now take a general element from Spin(3) and act on all three unit basis vectors e l7 e 2 ,e 3 . 
We can easily verify that the new elements e 11; e 2 2, e 33 provide another orthonormal basis with the 
same orientation: 

> ell:=rot(el,gSpin) ;e22 :=rot (e2,gSpin) ;e33 :=rot (e3,gSpin) ; 

ell := 2 (eps kc2 + c3 cl)e3 - (2 c2 2 + 2 c3 2 - 1) el + 2 {eps kc3 - c2 cl) e2 

e22 :=2(epsKcl - c3 c2) e3 - 2 {eps k c3 + c2 cl) el - (2 cl 2 -1 + 2 c3 2 ) e2 

e33 := -(2 c2 2 + 2 cl 2 - 1) e3 - 2 (-c3 cl + eps k c2) el - 2 (eps n cl + c3 c2) e2 

> el &c e2 + e2 &c el, el &c e3 + e3 &c el, e2 &c e3 + e3 &c e2; 

0, 0, 

> ell &c e22 + e22 &c ell, ell &c e33 + e33 &c ell, e22 &c e33 + e33 &c e22; 

0, 0, 

> el &w e2 few e3,ell few e22 &w e33; 

el 23, el 23 
Length of a vector v under the action Spin(3) is of course preserved: 

> Enorm(v) ,Enorm(Spin_action(v,gSpin) ) ; 

c3 2 + cl 2 + c2 2 , c3 2 + cl 2 + c2 2 

Example 1: Rotations in coordinate planes 

Let's define unit quaternions responsible for the rotations in the coordinate planes. These are 
counter-clockwise rotations when looking down the rotation axis. We will define a pure-quaternion 
basis consisting of {q^, q^, q^} in place of traditionally used {i, j, k}. 

> qi:=e23:qj :=el3:qk:=el2: 
ql2:=cos(alpha/2)*Id+sin(alpha/2)*'qk' ;#rotation in the xy-plane 

ql2 := cos(— a) Id + sin(— a) qk 

> ql3: =cos (beta/2) *Id+sin(beta/2)* 'qj ' ;#rotation in the xz-plane 

ql3 := cos(- /3) Id + sin(- ff) qj 

> q23:=cos(gamma/2)*Id+sin(gamma/2)*'qi ' ;#rotation in the yz-plane 

q23 := cos(- 7) Id + sin(- 7) qi 

Notice that to rotate by an angle na it is enough to find the n-th Clifford power of the appropriate 
quaternion and then apply it to a vector. 
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> ql2 &c ql2; #rotation by the angle 2*alpha 

cos(a) Id + el2 sin(a) 

> ql2 &c ql2 &c ql2; #rotation by the angle 3*alpha 

3 3 

cos(- a) Id + sin(- a) el2 

> ql2 &c ql2 &c ql2 &c ql2; #rotation by the angle 4*alpha 

cos(2 a) Id + sin(2 a) el2 
Let's see now how these basis rotations in the coordinate planes act on an arbitrary vector v 
aei + be 2 + ce 3 : 

> v:=a*el+b*e2+c*e3; 

V := a el + b e2 + c eS 
The norm of v is ||v|| = vv* and it can be defined in CLIFFORD as follows: 

> vlength: =sqrt (scalarpart (v &c star(v))); 



vlength := \/c 2 + a 2 + b 2 
Certainly, rotations do not change length. For example, let's rotate v by 913912 : 

> vl23:=rot(v,ql3 &c ql2) ; #rotation ql2 followed by ql3 

vl23 :— (— 6sin(/3) sin(a) + asin(/3) cos(a) + ccos(/3)) eS 

+ (— b sm(a) cos((3) — csin(/3) + ocos(/3) cos(a)) el + (fecos(a) + asin(a)) e^ 

> vlength: =sqrt (scalarpart (vl23 &c star (vl23) ) ) ; 



vlength := y/c 2 + a 2 + b 2 
Thus, the length of V123 is the same as the length of v. However, rotations do not commute. We 
will show that by applying quaternion 912913 to v and by comparing it with 9123 : 

> vl32:=rot(v,ql2 &c ql3) ; #rotation ql3 followed by ql2 

vl32 := (asin(/3) +ccos(/3)) e3 + (a cos(/3) cos(a) — csin(/3) cos(a) — 6sin(a)) el 
+ (6cos(a) — csin(/3) sin(a) + asin(a) cos(/3)) e2 

> clicollect(vl23-vl32) ; 

(— 6sin(/3) sin(a) + asin(/3) cos(a) — asin(/3)) e3 

+ (— 6sin(a) cos(/3) — csin(/3) + csin(/3) cos(a) + 6sin(a)) el 

+ (asin(a) + csin(/3) sin(a) — asin(a) cos(/3)) e^ 
As it can be seen, V123 ^ V132. 

Example 2: Counter-clockwise rotation by an angle a around the given axis 

In this example we will find a way to rotate a given vector v by an angle a in a plane orthogonal 
to the given axis vector axis = aiei + 0262 + 0363 vector. This rotation will be counter-clockwise 
when looking down the axis towards to origin (0, 0, 0) of the coordinate system. In order to derive 
symbolic formulas, we will assume that the symbolic vector axis has been normalized by defining 
A = ±\/l — a 2 — a 2 and axis = aiei + 0262 + Ae3. However, it won't be necessary for axis to be 
of unit length when its components are numeric. 

> alias (lambda=RootOf (-al~2-a2~2-_Z~2+l) ) : axis : =al*el+a2*e2+lambda*e3 ; 
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axis := al el + a2 e2 + A e# 

> simplif y(axis &c star(axis)); 

id 
Thus, in the symbolic case we will always have that ||axis|| = 1. Notice that in order to represent a 
rotation around axis we need to find a dual unit quaternion which we will call qaxis. It will need 
to be defined in such a way as to give a desired orientation for the rotation. Since we have opted 
for counter-clockwise rotations, we define qaxis = —axis e 123 where e 12 3 is a unit pseudo-scalar in 
C^o,3- In the following we will refer to the axis = axe! + a 2 e 2 + a 3 e 3 as the axis (ai, a 2 , 0,3). 

> qaxis :=axis &c (-el23) ; 

qaxis := A el2 + al e23 - a2 el3 
In Appendix 3 Reader can find a procedure qrot which finds the dual quaternion qaxis. The first 
three arguments to qrot are the (numeric or symbolic) components of the axis vector in the basis 
{ei,e 2 ,e3J while the fourth argument is the angle of rotation. For example, we can define various 
rotation quaternions: 

> ql00:=qrot(l,0,0,theta) ;#rotation about the axis (1,0,0) 

qlOO := cos(- 9) Id + sin(- 9) e23 

> q010:=qrot(0,l,0,theta) ;#rotation about the axis (0,1,0) 

qOlO := cos(- 9) Id - sin(- 9) el3 

> q001:=qrot(0,0,l,theta) ;#rotation about the axis (0,0,1) 

qOOl := cos(- 9) Id + sin(- 9) el2 

> ql01:=qrot(l,0,l,theta) ;#rotation about the axis (1,0,1) 

qlOl := cos(- 6) Id + sin(- 9) (-el2V2 + -V2 e23) 

> q011:=qrot(0,l,l,theta) ;#rotation about the axis (0,1,1) 

qOll := cos(- 9) Id + sin(- 9) (- el2 \/2 - - \/2 el3) 

> qll0:=qrot(l,l,0,theta) ;#rotation about the axis (1,1,0) 

qllO := cos(i 9) Id + sin(| 9) {)-^/2e23-)-^/2 el3) 

> qlll:=qrot(l,l,l,theta) ;#rotation about the axis (1,1,1) 

qlll := cos(- 9) Id + sin(- 9) (- el2 V3 - - V3 el3 + - VS e23) 

Zi Zi o o 

Let's try to rotate first a basis vector ei around various axes listed above by some angle a. 
In the next example, after we find a general formula for the components of the rotated el, we will 
verify it for various angles. 

> v:=el : vnew:=rot (v,ql00) ; #rotation around the (100) axis 

vnew := el 

> vnew:=rot (v,q010) ; #rotation around the (010) axis 

vnew :— — sin(6>) e3 + cos(8) el 
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> eval(subs(theta=Pi/2,vnew) ) ; 

-e3 

> vnew:=rot (v,q001) ; #rotation around the (001) axis 

vnew := cos(9) el + e2 sm(8) 

> eval(subs(theta=Pi/2,vnew) ) ; 

e2 

> vnew: =rot (v,ql01) ; #rotation around the (101) axis 

vnew := -- (-1 + cos(0)) e3 + - (cos(0) + 1) el + -V2e2sm{9) 

> eval (subs (theta=Pi/2, vnew) ) ; 

> vnew:=rot (v,q011) ; #rotation around the (Oil) axis 

vnew :=--\[2e3 sin(6») + cos(#) el + - y/2 e2 sin(6>) 

> eval (subs (theta=Pi/2, vnew) ) ; 

— e3V2+-e2V2 
2 2 



> vnew:=rot (v,qll0) ; #rotation around the (110) axis 

vnew := — — \[2 e3 sin 

> eval(subs(theta=Pi/2,vnew) ) ; 



vnew := -- \/2 e3 sin(6>) + - (cos(0) + 1) el - - (-1 + cos(9)) e2 



1^1 1 

~^eSV2+-el+-e2 

> vnew:=rot (v,qlll) ; #rotation around the (111) axis 

vnew := — (\/3sin(0) - 1 + cos(0)) e5 + - (2cos(6>) + 1) el 

+ - (V3sin(6>) + 1 - cos(0)) e2 

> eval (subs (theta=Pi/2, vnew)) ; 



-i(V3-l) e 5 + i e i+i(V3 + l)ei? 



> eval (subs (theta=Pi , vnew) ) 



2 1 2 

-eS--el + -e2 



Example 3: A general formula 

Finally, we derive a general formula for a rotation of an arbitrary vector v = t>iei + i>2e2 + 1*363 
around an arbitrary axis = aiei +0262 +0363 and by an arbitrary angle a. In the purely symbolic 
case we assume that the axis is of unit length, that is, 03 = A. 
> v:=vl*el+v2*e2+v3*e3; #an arbitrary vector 

v := vl el + v2 e2 + v3 e3 
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> qnew: =clicollect (rot (v,qrot (al ,a2, lambda, theta) ) ) ;# new vector after rotation 

qnew := (-A vl al cos(9) - Xv2 a2 cos(9) + al 2 v3 cos(8) - vl a2 sin(0) + A vl al 
+ Xv2 a2 + a2 2 v3 cos(0) + v2 al sin(0) + v3 - a2 2 v3 - al 2 v3)e3 + (vl cos(0) 

- al v2 a2 cos(0) - Xv3 al cos(9) + Xv3 al + al 2 vl - al 2 vl cos(0) 

- v2 Asin(0) + v3 a2sm(9) + al v2 a2)el + (a2 2 v2 + al vl a2 - v3 al sin(0) 

- a2 2 v2 cos(0) + A v3 a2 + vl A sin(0) - A v3 a2 cos(0) + v2 cos(0) 

- al vl a2cos(9))e2 

For example, let's rotate vector v = ei + 2e2 + 3e3 around the axis (1,2,3) by any angle a. 
Certainly, since the vector v is on the axis of rotation, it should not change: 

> clicollect(rot(el+2*e2+3*e3,qrot(l,2,3,alpha))) ; 

el + 2 e2 + 3 e3 
Let's rotate v = ei — 2e2 + 4e3 around the axis (2, —3, 4) by an angle a — 7r/4. 

> clicollect(rot(el-2*e2+4*e3,qrot(2,-3,4,Pi/4))) ; 

, 1 / — - 10 r- 96 x , 2 ,— ,- 48 19 _, 

Thus, in this section we have shown how easy it is to derive vector rotation formulas from vector 
analysis using elements of Spin(3) considered as unit quaternions. It has been very helpful to be 
able to embed Spin(3) in C£q,3- 



5.3 Degenerate Clifford algebra and the proper rigid motions SE(3) 

In this final section we will use the ability of CLIFFORD to perform computations in Clifford 
algebras of an arbitrary quadratic form including, of course, degenerate forms. We will consider 
the semi-direct product Spin(3) x M 3 that double covers the group of proper rigid motions SE(3). 



B 



We will follow the notation used in 27 , page 156, except that our basis vector that squares to 
will be e4 and not e. We begin by defining B as a degenerate diagonal form diag(— 1, — 1, — 1,0) 
of signature (0,3,1). Recall from the previous section that procedure star gives conjugation in 
C4,3,i = C(0,3,l). 

> dim:=4:n:=dim-l:eval(makealiases(dim)) :B:=diag(-l$n,0) ; #Selig's C(0,3,l) 

-10 
0-1 
0-10 

Let the vector basis in M 3 be stored in vbasis and let t an arbitrary vector in R 3 . 

> vbasis : =cbasis(n, 1) : t :=add(t . i*vbasis [i] , i=l . . n) ; 

t := tl el +t2e2 + t3 e3 
Elements of the form 1 + 1 e^ are invertible in C£o,3,i since u e^ and e^ u are nilpotent for any u 
in C£o,3,i- That is, the Jacobson radical J in C£o,3.i is generated by e^. 

> cinv(l+t &c e4) ;#symbolic inverse of 1+t & c e4 
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Id - tl e!4 - t2 e24 - t3 e34 



We will verify now statements made on page 156. Let k = \/ —c\ — c\ — c\ + 1 and e = ±1 as 
before. Then the most general element g in Spin(3) has the form: 

> alias (kappa=sqrt(-cl~2-c2~2-c3~2+l)) : alias (eps=RootOf (_Z~2-1) ) : 
gSpin:=eps*kappa*Id+c3*el2+c2*el3+cl*e23; 

gSpin := eps n Id + c3 el2 + c2 el3 + cl e23 
We consider a subgroup G of the group of units of C£o,3,i of the form g + ^tge 4 where g belongs 
to Spin(3) and t is a 1 -vector. Elements in G will be given by a procedure ge. 

> ge:=proc(g,t) RETURN(clicollect(simplif y(g+l/2*t &c g &c e4))) end: 
For the most general g in Spin(3) and tgR 3 , procedure ge gives: 

> 'ge(gSpin,t) '=ge(gSpin,t) ; 

ge(gSpin, t) = (-t2n eps - - tl c3 + -t3 cl) e24 +(-t3n eps - -t2 cl - -tl c2) e34 
+ (- t2 c3 + - t3 c2 + - tl k eps) el4 + (-- t2 c2 + - t3 c3 + - tl cl) el234 

£j Lj Zi £t £i Zi 

+ eps kH+c3 el2 + c2 el3 + cl e23 
First, let's verify Selig's statement that 

(g-^tge 4 )* = (g* + ig*te 4 ). (35) 

Notice that the left-hand-side in ( p5[ ) is just the conjugation of ge(gSpin, —t) while the right-hand- 
side is equal to ge(gSpin* 1 —t) where gSpin* denotes the conjugate of gSpin. 

> L: =clicollect (star (ge(gSpin,-t) ) ) : 
R:=simplify(star(gSpin)+l/2*star(gSpin) &c t fee e4) : 
simplif y(L-R) ; 


Next, we define the action of the group G on the subspace of C^o,3,i consisting of the elements of 
the form 1 + xe 4 as follows: 

(g+itge 4 )(l + xe 4 )(g-itge 4 )* = l + (gxg*+t)e 4 (36) 

where x, t € M 3 and g £ Spin(3). The above identity can be shown as follows. We define a 
procedure rigid which will give this action on R 3 : xntgxg'+t. 

> rigid: =proc(x,g,t) local p; 

if not evalb(x=vectorpart (x, 1) ) or not evalb(t=vectorpart (t , 1) ) then 

ERRORCx and t must be vectors') 
fi: 

if not type(g,evenelement) then ERRORCg must be even') fi; 
RETURN (clicollect (simplif y (cmul (g , x , star (g) ) +t ) ) ) 
end: 

This action will give us the rigid motion on IR 3 . Thus, we can compute the right-hand-side of (pO) 
by using rigid while the left-hand-side will be computed directly. 

> x: =add(x. i*vbasis [i] , i=l . .n) ; 

x := xl el + x2 e2 + x3 e3 
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> LHS:=simplify(ge(gSpin,t) &c (1 + x &c e4) &c star (ge(gSpin.-t) ) ) : 
RHS:=simplify(Id+rigid(x,gSpin,t) &c e4) : 
simplify(LHS-RHS); 


Finally, we will verify directly that the action xn gxg*+t is a rigid motion. If we denote by 
x p , y p the images of x, y under this action, we will need to show that ||x p — y p || = ||x — y|| in the 
Euclidean norm. We can compute the norm ||x — y|| by taking \/(x — y)(x — y)* in C^o,3,i, or by 
using a procedure distance. 

> y :=add(y . i*vbasis [i] , i=l . .n) : 

distance :=proc(x,y) sqrt (simplify (scalarpart (cmul(x-y , star (x-y) )) ) ) end: 
xp:=rigid(x,gSpin,t) :yp:=rigid(y ,gSpin,t) : 
evalb(distance(x,y)=distance(xp,yp) ) ; 

true 
Thus the action of G defined as xngxg*+t is a rigid motion in K 3 . In view of the presence of 
the radical in C^o,3,i, this group G is in fact a semi-direct product of Spin(3) and R 3 , that is of 
rotations and translations. It is well known of course that Spin(3) x M 3 doubly covers the group 
of proper rigid motions SE(3). We leave it as an exercise for the Reader to check in CLIFFORD 
that the composition of two rigid motions is a rigid motions. 

6 Summary 

The main purpose of this paper has been to show a variety of computational problems that can 
be approached with the symbolic package CLIFFORD. It is through an extensive experimentation 
with the package that the results reported in |7j were found; we have seen some of the computations 
that have led to them in Section |3j. In view of their complexity it is unlikely that they could have 
been performed by hand. Relation between the Clifford products in C£(g) and C£(B) through the 
Helmstetter's formula appears more clear once it has been checked with CLIFFORD. The SVD 
of a matrix as performed in Section [| is clearly feasible in the Clifford algebra language; however 
it is not clear if the approach presented there is the best. More study would need to be done 
here in order to possibly simplify the computations, perform them uniquely in the Clifford algebra 
language, and possibly better utilize the nilpotent-idempotent basis in the Clifford algebra rather 



than the Grassmann basis. In robotics applications presented in Section 5.3, CLIFFORD appears 
to be a very convenient tool to carry out practical computations in the low dimensional algebras 
such as C^o,3 and C£\^,i- 
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8 Appendix 1 

In addition to the main package CLIFFORD, in Section we have used new procedures from a 
supplementary package SUPPL. These procedures are: cinvg, cmulg, LCg, makeF, RCg, revg, and 
splitB. 

• Procedure cinvg finds the Clifford inverse with respect to the symmetric part g of the bilinear 
form B, that is, it finds the inverse (if it exists) of u in C£(g). 

• Procedure cmulg performs Clifford multiplication with respect to the symmetric part g of the 
bilinear form B : It gives the Clifford product uv for any two elements u and v in C£(g). 

9 

• Procedures LCg and RCg give the left and right contraction in C£(g). 

• Procedure makeF computes element F E f\ V defined in (0). It uses additional procedures 
pairs and mysign from the package SUPPL. 

• Procedure revg performs the reversion anti-automorphism " in C£(g). 

• Procedure splitB splits a bilinear form B in V into its symmetric part g and its antisym- 
metric part A which are returned as a sequence g,A. If B is purely symbolic and a second 
optional parameter (of any type) is used, in addition to g and A the output sequence contains 
two lists of symbolic substitutions that relate entries of g and A to the entries of B. If B is 
not assigned, in order for splitB to work it internally assigns a blank 9x9 matrix to B and 
then it calls itself. 

In Section we have used the following additional procedures from SUPPL: 

• Procedure cliexpand expands the given Clifford number u € C£(B) from the default Grass- 
mann basis to a Clifford basis consisting of un-evaluated Clifford products of the generators 
{ei, . . . , e„}. Procedure clieval converts from the Clifford basis back to the Grassmann basis. 

• Procedure reversion gives the reversion in C£(B). 

• Procedure clisolve2 solves equations of the type u = for the unknown parameters in u. 
Its code is shown in Appendix 3. 

• Procedure def B is needed to define a bilinear form B in an even-dimensional vector space V. 

• Procedure bexpand expands any element in the Hecke algebra Hf(n, q) in terms of the Hcckc 
& -basis. 

• Procedure alpha2 provides a Hecke algebra automorphism. 

In Section we have used the following additional procedures from a supplementary package AVSD. 

• Procedure phi provides an isomorphism between a matrix algebra and a Clifford algebra. 

• Procedure radsimplif y simplifies radical expressions in matrices and vectors. 
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• Procedure assignL is needed to write output from a Maple procedure eigenvects in a suit- 
able form, it sorts eigenvectors according to the corresponding eigenvalues, and it uses the 
Gram-Schmidt orthogonalization process, if necessary, to return a complete list of orthogonal 
eigenvectors. 

• Procedure climinpoly belongs to the main package CLIFFORD. It computes a minimal 
polynomial of any element of a Clifford algebra. 

• Procedure makediag makes a "diagonal" S matrix consisting of singular values. 

• Procedure embed embeds the given non-square matrix or a matrix of smaller dimensions into 
a 2 fe x 2 fc matrix of smallest k such that it can be mapped into a Clifford algebra. 

9 Appendix 2 

We display a Maple code of some more important procedures used in Section ||. 

Procedure bexpand expands any element in the Hecke algebra in terms of the Hecke b -basis: 

> suppl [bexpand] : =proc (X) local a, setbs, setbe, i,eq,T,sys,sol; 
option remember; 

setbs : = [Id, 'bl \ 'b2' , 'bl2' , 'b21 ' ,'bl21'] : 

setbe:=[Id,bl,b2,cmul(bl,b2),b21,cmul(bl,b21)] : 

eq:=CS(X-add(a[i]*setbe[i] ,i=l. .6)) ; 

T:=cliterms(eq) ; 

sys :={coef f s(eq,T)} ; 

sol :=solve(sys,{seq(a[i] ,i=l. .6)}) ; 

subs (sol, add(a[i] *setbs [i] , i=l . .6) ) ; 

end: 

Procedure alpha2 gives an automorphism in the Hecke algebra: 

> suppl [alpha2] :=proc(X) local a, setbb, setbe, setba, i,eq,T,sys,sol; 
option remember; 

setbe : = [Id,bl ,b2 , crnuKbl ,b2) ,b21 , crnuKbl ,b21)] : 

setba: = [Id,-l/q*reversion(bl) , -l/q*reversion(b2) , l/q~2*reversion(setbe [4] ) , 

l/q~2*reversion (setbe [5] ) , -l/q~3*reversion( setbe [6] )] : 

eq:=CS(X-add(a[i]*setbe[i] ,i=l. .6)) ; 

T:=cliterms(eq) ; 

sys :={coef f s(eq,T)} ; 

sol :=solve(sys,{seq(a[i] ,i=l. .6)}) ; 

subs (sol, add (a [i] *setba[i] , i=l . .6) ) ; 

end: 



Aliases needed in Section (3.3 



> alias (Hl=h [1] *q~3+2*h [6] *q-h [6] +h [5] *q~3-h [3] *q~2-2*h [5] *q~2-2*h [6] *q~2 

+h [5] *q+h [6] *q~3+h [4] *q~3+h [2] *q~3-fi [2] *q~2+h [4] *q+h [3] *q~3-2*h [4] *q~2) 
alias (H2=-2*h [6] *q+h [6] -h [4] *q+h [6] *q~2+h [2] *q~2+h [4] *q~2+h [5] *q~2-h [5] *q) : 
alias (H3=h [3] *q~2-h [5] *q-h [4] *q-2*h [6] *q+h [4] *q~2+h [6] *q~2+h [5] *q~2+h [6] ) : 
alias (H4=-h [6] +h [5] *q+h [6] *q) : 
alias (H5=h [6] *q+h [4] *q-h [6] ) : 
alias(H6=h[6]) : 



Aliases needed in Section ( |3.3|) 
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> alias (wl=q~4*K [4] ~2+3*q~3*K [4] "2-K [4] *q~3+4*q~2*K [4] ~2-K [4] *q~2 

+3*q*K [4] "2-K [4] *q-q+K [4] ~2-K [4] ) : 
alias (w2=K [4] *q~3+2*K [4] *q~2-q~2+2*K [4] *q+K [4] ) : 
alias (w3=-K [4] ~2+K [4] +K [2] +K [4] *K [5] *q-q~3*K [4] *K [2] -q*K [2] *K [4] 

+K [6] *q~2*K [4] +K [4] *q~4*K [6] +K [6] *q~3*K [4] +K [4] *K [6] *q-q~2*K [2] *K [4] 

+q~2*K [4] *K [5] -K [2] *q-K [5] *q+K [4] *q~2-K [4] *q+q~4*K [4] ~2-q~2*K [4] ~2 

-q*K [4] "2-K [6] *q+K [6] *q~2-K [4] *K [2] ) ; 
alias (w4=K [6] *q~3+K [4] *q~3-K [6] *q~2-K [5] *q~2-K [6] *q-K [5] *q+K [6] +K [4] ) : 
alias (w5=K [4] *q+K [2] *q-K [6] -K [5] ) : 

alias (w6=K [6] *q~2+K [4] *q~2-2*K [6] *q-K [5] *q-K [4] *q+K [6] +K [4] ) : 
alias (w7=K [4] *q~3+2*K [4] *q~2+q+2*K [4] *q+K [4] ) : 
alias (w8=q~5*K [4] ~2+3*q~4*K [4] "2+4*q~3*K [4] ~2+K [4] *q"3+3*q"2*K [4] ~2 

+K [4] *q~2+q*K [4] "2+K [4] *q-l+K [4] ) : 
alias (w9=K [4] *q~3+K [2] *q'3-2*K [6] *q~2-K [5] *q"2-K [4] *q"2+K [5] *q 

+2*K [6] *q+K [4] *q-K [4] -K [6] ) : 

Aliases needed in Section ( |5.1[ ): 

> alias (kappal=RootDf (_Z~2-l+x3"2+x2"2+x4"2) ) : 
alias (kappa2=RootDf (_Z~2-l+x5~2+x3~2+x2~2) ) : 

alias (kappa3=Root0f (-x6~2+x5~2*x6~2+x4~2*x5~2+x6~4+x2~2*x6~2+x4~2*x6~2+_Z~2) ) : 
alias (kappa4=Root0f(-x7"2+x5"2*x7"2+x3"2*x7"2+x6"2*x7"2+x7"4+x4"2*x5"2 

-2*x4*x5*x6*x3+x6~2*x3~2+x4~2*x7~2+_Z~2) ) : 
alias (kappa5=Root0f ( (x7~2+x8~2) *_Z"2+ (2*x5*x4*x7-2*x3*x6*x7) *_Z-x8~2+ 

+X6"2*x3"2+x8"4+x4~2*x5~2-2*x4*x5*x6*x3+x6~2*x8~2+x7~2*x8~2+x5~2*x8~2 

+x3"2*x8"2+x4"2*x8"2)) : 
alias (eps=RootOf (_Z"2-1)) : 
alias (lambdal=RootOf (_Z~2+x2~2+x3~2-l) ) : 
alias (lambda2=Root0f (_Z"2-l+x3"2) ) : 
alias (lambda3=Root0f (_Z"2-l+x5"2) ) : 

alias (lambda4=Root0f (-x6"2+x5"2*x6"2+x4"2*x5"2+x6"4+x2"2*x6"2+x4~2*x6"2+_Z~2) ) : 
alias (lambda5=Root0f (_Z"2-x6"2+x5"2*x6"2+x6"4) ) : 
alias (lambda6=Root0f (-x7~2+x5~2*x7~2+x3~2*x7~2+x6~2*x7~2+x7~4+x4~2*x5~2 

-2*x4*x5*x6*x3+x6~2*x3~2+x4~2*x7~2+_Z~2) ) : 
alias (lambda7=Root0f (_Z"2-x7"2+x5"2*x7"2+x6"2*x7"2+x7"4) ) : 
alias (lambda8=Root0f ( (x7"2+x8"2) *_Z~2+ (2*x5*x4*x7-2*x3*x6*x7) *_Z-x8"2 

+X6"2*x3"2+x8"4+x4~2*x5~2-2*x4*x5*x6*x3+x6~2*x8~2+x7~2*x8~2+x5~2*x8~2 

+x3"2*x8"2+x4"2*x8~2)) : 
alias (lambda9=Root0f (_Z"2*x8"2-x8"2+x8"4+x3"2*x8"2+x4"2*x8"2) ) : 

10 Appendix 3 



We display code for two procedures rot and qrot that were needed in Section 5.2. Procedure rot 
performs a rotation of a vector by a quaternion through a certain angle. 
> rot :=proc(v,q) local qs ; 

qs :=star (q) : 

RETURN (map (factor , clicollect (simplify (cmul (q, v , qs) ) ) ) ) 

end: 

Procedure qrot finds a unit quaternion that is dual to the rotation axis vector axis. 
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> qrot : =proc(pl ,p2,p3,theta) local bas, c,e,k,l, i,q,n;global qaxis; 
if type(pl ,name) or type (p2, name) then 

RETURN (cos (theta/2) *Id+sin(theta/2) *qaxis) f i ; 
if evalb(simplify(pl~2+p2~2+p3~2)=l) then 

q:=simplif y (subs({al=pl,a2=p2} , qaxis)) f i; 
n:=sqrt(pl~2+p2~2+p3~2) : 
if n=0 then ERRORCaxis vector must be a non-zero vector') elif 

has(n,RootOf ) then n: =max(allvalues(n) ) fi; 
q:=simplif y( subs ({al=pl/n,a2=p2/n}, qaxis) ) ; 
bas: = [Id,el2,el3,e23] :c: = [] :k:=0: 
for i from 1 to 4 do 1 :=coef f (q.bas [i] ) ; 

if has(l,RootOf ) then k:=i fi: 
c:=[op(c) ,1] 
od: 
if k<>0 then e :=allvalues (c [k] ) ; 

if p3>0 then c : =subsop(k=max(e) , c) elif 

p3<0 then c : =subsop(k=min(e) , c) else ERR0R( 'p3=0' ) fi; 

fi; 

q:=add(c[i]*bas[i] ,i=l. .4) ; 

RETURN (cos (theta/2) *Id+sin (theta/2) *q) 

end: 

Procedure clisolve2 was used extensively throughout this paper to solve linear equations in a 
Clifford algebra C£(B). 

> suppl [clisolve2] :=proc(eq, indet) local i,T, vars, sol, sys; 
if type (indet , list) then vars :=convert (indet , set) 

else vars : =select (type , indets (indet) , indexed) 

fi; 

T:=cliterms(eq) ; 

sys :={coeffs(clicollect (simplify (eq) ) ,T)} ; 

sol:= [solve(sys , vars)] ; 

if type (indet, list) then RETURN (sol) 

else RETURN ( [seq(subs(sol [i] , indet) , i=l . .nops(sol) )] ) ; 

fi; 

end: 
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